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I .   INTRODUCTION 

A.   STATEMENT  OF  THE  PROBLEM 

In  the  past,  numerous  highway  bridges  have  suffered  extensive 
damages  due  to  strong  motion  earthquake  [1]  .   The  older  bridges,  con- 
sisting of  single  or  multiple  simple  truss  or  girder  spans  supported  on 
massive  piers  and  abutments,  were  particularly  vulnerable  to  the  action 
of  strong  ground  motions.   Seismic  damages  were  most  commonly  caused  by 
foundation  failures  resulting  from  excessive  ground  deformation  and/or 
loss  of  stability  and  bearing  capacity  of  the  foundation  soils.   As  a 
direct  result,  the  substructures  often  tilted,  settled,  slid,  or  even 
overturned;  thus  severe  cracking  or  complete  failure  was  often 
experienced.   These  large  support  displacements  also  caused  relative 
shifting  of  and  damage  to  the  superstructures,  induced  failures  within 
the  bearing  supports,  and  even  caused  spans  to  fall  off  their  supports. 
It  is  significant  to  note  that  very  little  damage  occurred  to  these 
older  structures  as  a  direct  result  of  structural  vibration  effects. 

Certain  types  of  modern  highway  bridges  may,  on  the  other  hand,  be 
quite  susceptible  to  damage  from  strong  ground  vibration  effects.   This 
fact  became  very  evident  during  the  San  Fernando  earthquake  of  February 
9,  1971,  when  numerous  reinforced  concrete  highway  bridges  which  were 
designed  by  the  traditional  elastic  approach;  i.e.,  by  the  equivalent 
static  seismic  coefficient  method  similar  to  that  previously  adopted  for 
buildings,  suffered  severe  damages.   Because  of  this  experience,  it  was 
immediately  apparent  that  the  seismic  requirements  used  were  inadequate 
and  that  the  design  methodology  should  be  critically  examined.   Action  was 
quickly  taken  following  the  earthquake  to  correct  certain  design 
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deficiencies;  however,  the  basic  static  approach  to  design  still  remains 
in  effect. 

Recognizing  the  urgent  need  for  both  theoretical  and  experimental 
research  related  directly  to  seismic  effects  on  bridge  structures ,  an 
investigation  entitled  "An  Investigation  Of  The  Effectiveness  of  Existing 
Bridge  Design  Methodology  in  Providing  Adequate  Structural  Resistance  to 
Seismic  Disturbances"  was  initiated  in  1971  within  the  Earthquake 
Engineering  Research  Center,  University  of  California,  Berkeley,  under 
the  sponsorship  of  the  U.S.  Department  of  Transportation,  Federal  Highway 
Administration.   This  investigation  consists  of  the  following  six  phases: 

(1)  A  thorough  review  of  the  world's  literature  on  seismic  effects 
on  highway  bridge  structure  including  damages  to  bridges  during 
the  San  Fernando  earthquake  of  February  9,  1971. 

(2)  An  analytical  investigation  of  the  dynamic  response  of  long, 
multiple- span,  highway  overcrossings  of  the  type  which  suffered 
heavy  damages  during  the  1971  San  Fernando  earthquake. 

(3)  An  analytical  investigation  of  the  dynamic  response  of  short, 
single  and  multiple  span  highway  overcrossings  of  the  type 
which  suffered  heavy  damages  during  the  1971  San  Fernando 
earthquake . 

(4)  Detailed  model  experiments  on  a  shaking  table  to  provide 
dynamic  response  data  similar  to  prototype  behavior  which  can 
be  used  to  verify  the  validity  of  theoretical  response 
predictions. 

(5)  Correlation  of  dynamic  response  data  obtained  from  shaking 
table  experiments  with  theoretical  response  and  modification 
of  analytical  procedures  as  found  necessary. 


(6)   Preparation  of  recommendations  for  changes  in  seismic  design 
specifications  and  methodology  as  necessary  to  provide  ade- 
quate protection  of  reinforced  concrete  highway  bridges 
against  future  earthquakes. 

Final  reports  coverning  Phases  1,  2,  and  3  have  been  published  [1,2,6] 
and  the  final  report  coverning  Phase  4  is  now  nearing  completion.   The 
present  final  report  covers  Phase  5  and  the  final  report  on  Phase  6  will 
be  completed  during  summer  and  fall  of  1977. 

B.   OBJECTIVES  AND  SCOPE 

The  primary  objectives  of  the  investigation  in  Phase  5  are  to  cor- 
relate the  dynamic  responses  of  the  experimental  model  measured  during 
shaking  table  tests  conducted  in  Phase  4  with  the  corresponding  responses 
determined  using  the  analytical  model  and  numerical  procedures  developed 
in  Phase  2. 

To  achieve  these  objectives,  the  analytical  model  and  numerical 
procedures  defined  and  developed  in  Phase  2  have  been  re-examined  for 
use  in  predicting  the  seismic  response  of  the  experimental  model  structure, 
Modifications  of  the  analytical  model  simulating  nonlinear  dynamic 
behavior  of  expansion  joints  and  the  analytical  procedure  for  numerical 
integration  of  equations  of  motion  have  been  found  necessary.   These 
changes  have  been  incorporated  into  the  nonlinear  response  analysis 
computer  program  NEABS. 

Experimental  test  results  obtained  in  Phase  4  were  reviewed  to 
select  suitable  excitation  tests  for  correlation  purposes.   Since  it  was 
found  that  the  model  structure  is  more  susceptible  to  damage  when  the 
excitation  is  oriented  in  the  transverse  direction,  test  results 
obtained  under  transverse  and  under  simultaneous  transverse  and  vertical 
excitations  were  used  for  correlation  purposes. 
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Correlation  of  the  measured  experimental  displacement  response  with 
the  corresponding  response  predicted  by  linear  and  nonlinear  analyses 
has  been  performed  for  three  types  of  seismic  excitation  (Series  IV 
excitation  tests),  namely,   (1)  low  intensity  excitation  in  the  trans- 
verse direction,   (2)  high  intensity  excitation  in  the  transverse  direc- 
tion, and   (3)  high  intensity  excitation  in  both  the  transverse  and 
vertical  directions.   The  linear  response  analysis  computer  program  BSAP 
and  the  nonlinear  response  analysis  computer  program  NEABS  as  revised 
have  been  used  for  these  correlations. 

Finally,  a  parameter  study  of  the  dynamic  response  was  conducted 
using  the  nonlinear  analytical  model  and  analysis  procedures. 

Chapter  II  of  this  report  discribes  the  modifications  introduced 
into  the  nonlinear  analytical  model  of  expansion  joints  which  allow  pro- 
per simulations  of  the  impact  phenomenon,  slippage  with  Coulomb  type 
friction,  and  separation  under  the  action  of  elastoplastic  joint  res- 
trainer  tie  bars  acting  in  tention  only.   Chapter  III  describes  certain 
modification  to  the  analytical  procedures  including  an  equilibrium  cor- 
rection introduced  through  iteration  and  subdivision  of  time  intervals 
into  a  specified  number  of  equal  sub- intervals.   Chapter  IV  describes 
the  experimental  model  and  test  results  which  were  used  to  formulate  the 
analytical  model  used  in  the  correlation  study.   Chapter  V  describes  the 
procedure  used  to  formulate  a  basic  analytical  model  which  reproduces 
the  observed  dynamic  characteristics  of  the  experimental  model  under  the 
initial  low  intensity  excitations.   Chapter  VI  presents  the  correlation 
results  for  the  three  seismic  excitations  having  various  intensities  of 
table  motion.   Chapter  VII  gives  a  discussion  of  the  correlation  results 
and  Chapter  VIII  summarizes  the  conclusions  drawn  from  the  investigation. 


Finally,  numerical  results  obtained  from  the  parameter  studies  of  impact 
and  Coulomb  friction  are  presented  in  order  to  show  the  applicability  of 
the  analytical  model  and  numerical  procedures  defined  in  Chapters  II  and 
III. 


II.   ANALYTICAL  MODEL  OF  EXPANSION  JOINTS 

A  nonlinear  analytical  model  for  simulating  the  dynamic  behavior  of 
expansion  joints  was  formulated  by  Tseng  and  Penzien  [2,4],   It  includes 
relative  translational  and  rotational  degrees  of  freedom,  elastoplastic 
joint  restrainer  tie  bars  acting  in  tension,  impact,  and  Coulomb  type 
friction  with  slippage.   Although  the  original  analytical  model  of 
expansion  joints  was  basically  correct,  certain  modifications  were  needed 
to  make  it  fully  reproduce  the  characteristics  of  the  model  bridge,  struc- 
ture.  The  modifications  introduced  into  the  original  analytical  model 
and  into  the  nonlinear  analysis  computer  program  NEABS  [2]  were  the 
following: 

(1)  The  impact  spring  was  modified  in  such  a  way  that  it  approx- 
imates elastic  collisions. 

(2)  The  Coulomb  type  friction  force  was  approximated  by  an  elasto- 
plastic hysteretic  model. 

The  subsequent  discussion  in  this  chapter  concentrates  on  the  above 
modifications  to  the  analytical  model  of  the  expansion  joints  with  some 
mention  of  the  basic  equations  needed  to  characterize  their  nonlinear 
behavior. 

A.   PIECEWISE  LINEAR  EXPANSION  JOINT  STIFFNESS  [2,4] 

The  force  vs.  displacement  relation  of  the  idealized  expansion  joint 
as  shown  in  Fig.  2.1  can  be  characterized  using  an  expansion  joint 
coordinate  system  as  defined  in  Fig.  2.2,  i.e., 
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in  which  \p     represents  the  skew  angle  of  the  deck. 
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and  the  corresponding  expansion  joint  force  F  defined  by 
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Thus,  the  idealized  expansion  joint  can  be  characterized  by  a  stiffness 
matrix  k  relating  Af  to  Au  ,   i.e., 
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The  stiffness  matrix  k    must  be  transformed  to  the  local  nodal 

coordinate  system  using  Eqs.  (6)  and  (7).   This  transformation  results 

in  the  relation 
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Matrix  k    is  the  stiffness  matrix  which  relates  the  vector  of  local 

nodal  force  increments  As  to  the  vector  of  local  nodal  displacement 

increments  Ar  at  time  t,   and  which,  after  local  to  global  coordinate 

transformation,  can  be  used  to  assemble  the  total  stiffness  matrix  K 

for  the  complete  structural  system. 

Stiffness  matrix  k  defined  in  Eq.  (11)  represents  contributions 

from  two  sources.   One  contribution  is  a  stiffness  matrix  k   of  the 

idealized  joint  without  collisions  and  another  contribution  is  matrix 

k   that  represents  the  effect  of  collisions,  i.e., 
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The  stiffness  coefficients  of  k   can  be  written  in  the  form 
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where  k   is  the  elastic  stiffness  of  the  spring  which  resists  the 
relative  transverse  displacement,   k   is  the  elastic  stiffness  of  the 
springs  which  resists  the  relative  vertical  displacement  at  points  A 
and  B ,   and 
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in  which  N   is  the  number  of  tie  bars  and  k  .   is  the  instantaneous 
T  Ti 

stiffness  of  the  ith  tie  bar  at  time  t. 

The  stiffness  coefficients  of  k   is  described  in  the  following 
paragraph . 

B.   COLLISIONS 

Longitudinal  collisions  are  defined  to  take  place  at  points  A  and 

B  when  the  relative  displacement  between  the  two  end  diaphrams  close 

the  joint  gap  A   with  a  nonzero  relative  velocity.   At  the  instant 
G 

collision  takes  place,  the  longitudinal  impact  springs  having  large 

stiffness  k   which  are  attached  to  one  end  diaphram  leaving  a  small 

gap  A   with  the  other  end  diaphram  start  to  resist  the  motion.   A 
G 

collision  is  completed  when  rebound  occurs  and  the  relative  displace- 
ment between  the  two  diaphrams  becomes  equal  to  the  joint  gap  A  . 

G 

The  contact  forces  acting  at  points  A  and  B  can  be  written  as 
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1     u   +  A„  <  0 
Ax    G 

0     u    +  A      >  0 
Ax    G  — 


1     u  +  A   <  0 
Bx    G 

0     u„  +  A„  >  0 
Bx     G  — 


(18) 


If  it  is  assumed  that  (u   +  A  )  and  (u   +  A  )  do  not  change  sign 

Ax     G         Bx     G 

during  a  time  interval  At,   the  change  of  contact  force  during  a  time 
interval  can  be  expressed  as 


Ap.t  =  kT  <  u   +  A^  >  Au 

AI       I  Ax    G     Ax 

AP„T  =   k  <  u        +  A  >  &u 

BI       I  Bx    G     Bx 


(19) 


Let  Af   represent  the  incremental  contact  force  vector  in  the 
relative  expansion  joint  coordinate  system  during  a  time  interval  At 
as  defined  by 


^ 


Ap. 


Af1  =1 


AI 
0 

0 


Ap. 


BI 
0 


0 

^-  J 


Then,   Af   can  be  expressed  as 


(20) 


Af1  =  k  Au 


(21) 


in  which 
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*2 


k     <   u        + 
I           Ax 

A 
G 

> 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

kI 

< 

UBx 

+ 

A 
G 

> 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

(22) 


C.   COULOMB  FRICTION  FORCES 

Coulomb  friction  forces  are  developed  at  contact  points  A  and  B 
when  the  expansion  joint  undergoes  longitudinal  relative  displacement 
and  when  the  vertical  contact  forces  are  compressive.   The  friction  at 
each  point  A  and  B  always  acts  in  the  direction  opposite  to  the 
relative  velocity  as  a  pair  of  self  equilibrating  forces.   When  the 
expansion  joint  does  not  undergo  longitudinal  relative  displacement, 
each  friction  force  can  have  a  magnitude  anywhere  between  its  maximum 
and  minimum  values.   The  magnitude  of  a  friction  force  depends  on  the 
magnitude  of  other  forces  acting  on  the  expansion  joint.   Such  char- 
acteristics can  be  represented  by  a  rigid-plastic  hysteretic  force- 
relative  displacement  model  as  shown  in  Fig.  2.3.   It  should  be  noted 
here  that  the  model  of  the  Coulomb  friction  force  represented  in 
Fig.  2.3  is  mathematically  the  same  as  the  model  of  a  rigid-plastic 
restoring  force.   Plastic  deformations  in  the  latter  model  correspond 
to  slippages  in  the  former  model.   In  order  to  avoid  numerical 
instability  caused  by  a  sudden  change  in  the  Coulomb  friction  force  at 
zero  relative  velocity,  the  force  is  modeled  by  the  elastoplastic 
hysteretic  force-relative  displacement  model  shown  in  Fig.  2.4.   If  the 
slope  of  OA  in  Fig.  2.4  is  taken  large  enough,  then  both  models  of 
Figs.  2.3  and  2.4  are  practically  the  same. 
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According  to  the  assumption  described,  the  Coulomb  friction  forces 
acting  at  points  A  and  B  at  time  t  can  be  expressed  as 


'Ax 


-{ 


-V  <  F,   >   F„ 
Az     Az 


kC  <  FA   >  (u   -  v?    ) 
Az     Ax    Ax 


^V<*te>  FAZ! 


for   u   <   u 

Ax  _  Ax 


u 


Ax 


S      E   ^      ^   S      E 
for   u   -u   <u   <u   +u 

Ax    Ax    Ax    Ax    Ax 


^   S     E 
for   u   >  u   +  u 

Ax  —  Ax    Ax 


(23) 


r 


-V  <  F    >   F 

Bz     Bz 


-»    =  (  k°  <  F„   >  (uD   -  uf  ) 

Bx     A       Bz     Bx    Bx 


Bz     Bz 


for   u   <  u   -  u 

Bx  —  Bx    Bx 


_      S     E   .      .   S     E 
for   u    -u    <u    <u    +u 

Bx    Bx    Bx    Bx    Bx 


for   u    >  u    +  u 

Bx  ~  Bx    Bx 


S         S 
where  v  is  a  constant  coefficient  of  Coulomb  friction,   u    and  u 

Ax        Bx 

E 
are  the  current  slippages  at  points  A  and  B,   respectively,  and  u 

AX 

E 
and  u    are  the  elastic  deformations  at  points  A  and  B,   respectively, 
Bx 

as  given  by 


u    =  v  If 

Ax 


Az 


A, 


u     =   V 
Bx 


F«  A„ 
Bz    C 


(24) 


where 


<   F         > 
Az 

=    ■ 

1 
0 

F,       <    0 
Az 

F,      >   0 

Az  — 

<   F         > 
Bz 

= 

1 
0 

F„      <   0 
Bz 

FD      >   0 
Bz  — 

(25) 


If  it  is  assumed  that  the  vertical  compressive  contact  forces  F    and 

S      E  S      E 

F_   are  constant  and  deformations  (u„   -  u_   +  u,  )  and  (u„  -  u_  +  un  ) 
Bz  Ax     Ax  -   Ax         Bx     Bx  -   Bx 
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do  not  change  signs  during  a  time  interval  At,   then  the  changes  of 
Coulomb  friction  force  during  the  time  interval  can  be  expressed  as 


r 


k   <  p    >  Au 

Az     Ax 


ACAX  "  <  ° 


v.. 

r 


k  <  F„  >  A» 

Bz      Bx 


Ac    =  (  0 


Bx 


S     E   .      .   S     E 
for    u    -  u_   <  u^   <  u    +  u, 
Ax     Ax     Ax     Ax     Ax 

S      E 
for    u,   <  u    -  u    or 

Ax  —  Ax    Ax 


.   S      E 
u    >  u    +  u 
Ax  ~~  Ax     Ax 


S      E  S      E 

f°r    UBx  ~  UBx  <  UBx  <  uBx  +  uBx 


S      E 
for    u„   <  u^   -  u^   or 
Bx  —  Bx    Bx 


(26) 


^ 


S      E 
'Bx  "  aBx   UBx 


,-C 


Let  AP   represent  the  incremental  Coulomb  friction  force  vector 


in  the  expansion  joint  coordinate  system  during  the  time  interval  At; 


,-C 


then,   AP   can  be  expressed  as 


AfC 


€ 


€ 


r  -\ 

-Ac, 

Ax 
0 


r 


€  ■< 


AC. 


Ax 
0 


-AC 


K  •< 


Bx 
0 


v_ 


J 


Ac 


Bx 


*\ 


\    (27) 


V_ 


J 


Thus  AP   can  be  related  to  the  vector  of  incremental  expansion  joint 

-C 
displacement  Ar  by  a  matrix  k   in  the  form 


—c  -c       - 

Apu  =  k   Ar 


(28) 


-C 
where  k   is  a  function  of  r  at  time  t. 
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-C 
From  Eqs.  (9)  and  (10) ,  k   can  be  written  in  the  form 


*t 


*3 


h 


-*3 


k-3 


(29) 


in  which  k   is  a  matrix  defined  in  the  relative  expansion  joint 
coordinate  system  which  relates  the  change  of  Coulomb  friction  force 
vector  during  a  time  interval  At  to  the  change  of  relative  expansion 
joint  displacement  Au.   From  Eqs.  (26)  and  (27),   k   can  be  written 
as 


where 


*3 


A 
0 

0 

0 

0 

0 


0 
0 
0 

o 


(30) 


r 


k   <  F,   > 
Az 


kC  -  <  0 
A     ^ 


*^~ 


r 


4  =  <° 


for 
for 


k  <  F„   >    for 
Bz 


for 


S     E  .  ,     S   ,   E 

Ax    Ax    Ax    Ax    Ax 

.   S      E 
u,   £  u,   -  u_   or 
Ax    Ax    Ax 

>   S   4.  E 
u,  >   u„   +  u_ 
Ax  —  Ax    Ax 

(31) 


S      E   .      ,      S      E 
u   -  u    <  u    <  u„  +  u_ 
Bx    Bx    Bx     Bx    Bx 

.   S      E 
UBx  -  UBx  "  UBx  °r 

> uS   + E 
Bx  —  Bx    Bx 
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-C 
Finally,  the  matrix  k   must  be  transformed  to  the  local  coordinate 

system  using  Eqs.  (6)  and  (7).   This  transformation  results  in  the 

relation 


C      C  C      T  -C 

AP   =  k   Ar    ,    k   =  A   k   A  (32) 

C 
Matrix  k   relates  the  vector  of  local  nodal  Coulomb  friction  force 

C 
increment  AP   to  the  vector  of  local  nodal  displacement  increment 

Ar.   After  local  to  global  coordinate  transformation,  it  can  be  used 

to  assemble  the  total  stiffness  matrix  K   together  with  the  stiffness 

,  EJ 
matrix  k    ,  ,.   ,  .   _    .,_. 
-t   defined  in  Eq.  (13)  . 
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eAM9, 


.   % 


LOCAL    NODAL    COORDINATE    SYSTEM    r 


EXPANSION    JOINT    COORDINATE    SYSTEM    r 


FIG.  2.2   COORDINATE  SYSTEM  FOR  EXPANSION  JOINT  [2] 
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^AX^BX^ 


Uax(Ubx) 


FIG.  2.3   RIGID-PLASTIC  HYSTERETIC  MODEL  FOR  COULOMB  FRICTION 


A  Fax  (^bx) 


Uax(Ubx) 


FIG.  2.4   ELASTOPLASTIC  HYSTERETIC  MODEL  FOR  COULOMB  FRICTION 
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III.   DYNAMIC  ANALYSIS  PROCEDURE 

A.   INCREMENTAL  EQUATION  OF  MOTION 

The  equation  of  motion  for  an  n  degree  of  freedom  system 
representing  dynamic  equilibrium  at  time  t  can  be  expressed  as 

It     +  *t  +  *t  =  h  (33) 

IDS 
where,   F  ,      F  ,   F   and  R   are  the  nodal  inertia  force  vector,  nodal 

damping  force  vector,  nodal  restoring  force  vector  and  nodal  external 

force  vector  at  time  t,   respectively.   When  the  structure  is  nonlinear, 

Eq.  (33)  is  written  for  time  t  +  At  in  the  form 

(?t  +  *£)*(£  +  ^t)+(^  +  ASt)=  strtt        (34) 

ID  S 

where,   AF  ,   AF   and  AF   are  the  changes  of  the  nodal  inertia 

forces ,  nodal  damping  forces  and  nodal  restoring  forces  during  a  time 

interval  At,   respectively.   These  changes  of  the  forces  are  assumed 

to  be  given  by 


AF   =  M  Au  ;  AF°  =  C  Au  ;  AFS  =   K  Au  (35) 


where,   M,   C  and  K   represent  the  constant  mass  matrix,  constant 
viscous  damping  matrix  and  tangent  stiffness  matrix  at  time  t,   respec- 
tively.   Vectors  Au  ,   Au   and  Au   are  the  changes  of  nodal 
accelerations,  nodal  velocities  and  nodal  displacements  during  time 
interval  At  as  defined  by 
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Au  =  u  -  u 
-t     -t+At     -t 

Aut  =  Ht+At  -  ut  (36) 

Au  =  u   .  -  u 
-t     -t+At     -t 


Substituting  Eq.  (35)  into  Eq.  (34)  ,  one  can  obtain  the  equation 
of  motion  in  incremental  form  as 


M  Au   +  C  Au.  +  K  Au.  =  R.  ...  -  M  u.  -  C  u.  -  pf  (37) 

-   -t   -   -t   -t  -t   -t+At   -  -t   -  -t   -t 


The  nodal  external  force  vector  R   .    is  assumed  to  be  derived 

-t+At 

from  two  types  of  loading  [2] ;  an  applied  dynamic  force  P   .   and  an 

inertia  force  due  to  ground  motion  u       ,    i.e., 

-g(t+At) 


*t+At  =  *t+At  "  *  5  Hg(t+At,  (38) 


in  which  B  is  a  matrix  of  ground  motion  influence  coefficients. 

The  tangent  stiffness  matrix  K   consists  of  time  independent 

coefficients  K   which  are  assembled  before  excitation  based  on  the 

N 
lnxtial  stress  stage  and  time-dependent  coefficients  K   which  change 

in  each  time  internal ,  i.e.. 


L      N 
K   =  K   +  K™  (39) 


where 


LEL  NEL 

j=l   "(3)      A   -(D) 


VN  _  NyL   N  N     _ 

-t     *   -t(j)  '    -t(j)  "  -t(j)  "  5(j) 


(40) 
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in  which  LEL  and  NEL  represent  numbers  of  linear  and  nonlinear 

elements  in  the  system.   Matrices  k . . .   and  k  , . ,   represent  time 

-(D)       -t(;j) 

independent  stiffness  coefficients  for  the  jth  linear  or  nonlinear 

element  and  time  dependent  stiffness  coefficients  for  the  jth  nonlinear 

S 
element,  respectively.   From  Eq.  (39) ,   F   can  be  written  as 

S      SL      SN      L         N?L    SN 

?t   -  *t  +  Et  =5  2t  +  lx    £t(j,        («) 

in  which  f   .    represents  time  dependent  nodal  restoring  forces  for 

the  jth  nonlinear  element.   It  is  desirable  to  calculate   f  ...  (j  £  NEL) 

-t(u)  J  v- 

directly  from  current  element  deformations  in  order  to  avoid  accumulating 
errors  during  the  numerical  computation  [16] . 

Ihe  damping  matrix  is  assumed  to  be  of  the  form 


C  =  a  M  +   B  KL  (42) 


showing  one  part  proportional  to  the  mass  distribution  and  the  other 
part  proportional  to  the  initial  stiffness  distribution  before  excitation. 

B.   STEP- BY-STEP  INTEGRATION 

Newmark1  s  generalized  acceleration  method  assumes  the  following 
approximations  for  the  nodal  velocities  and  displacements  [12] 


VAt  =  Ht  +  [d  -  ■«)'  St  +  «  5t+At]  At 

=  at  +  ut  At  +  [(±-o)   ut  +  aut+AJ  At 


(43) 


^t+At 


where  parameters  6  and  0     can  be  chosen  to  give  the  required 
integration  stability  and  accuracy.   When  6  =  1/2  and  0   =  1/6,   the 
approximations  correspond  to  the  linear  acceleration  method  and  when 
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<5  =  1/2  and  O   =  1/4,   they  correspond  to  the  constant  acceleration 
method. 

Using  Eq.  (36) ,  the  approximations  given  by  Eq.  (43)  can  be 
expressed  in  the  incremental  form 


Au   =cAu-Cu-Cu 
-t      1  -t    2  -t    3  -t 


Au   =  C.  Au  -  C_  u  -  r     ti 
-t      4  -t    5  -t    6  -t 


(44) 


where 


a  At' 

6 
aAt 


aAt  '   3 


C  = 
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=  (M 


1  At 


(45) 


Substituting  Eqs.  (41)  and  (44)  into  Eq.  (37)  ,  one  can  obtain 


C,  M  +  C.  C  +  K 
1  -    4  -   -t 


+  M 


Aut  =  ?t+At  -  m  ut  -  c  ut  -  kl  ut  -  rf 


2t  +  °3  5tJ  +  S  [C5  ?t  +  C6  »t 


(46) 


Using  Eq.  (42),  i.e.,   C  =  a  M  +  3  K  ,   and  introducing  the  following 
normalized  constants. 


c?  =  C  ■  +  a  c4, 


c8  =3c4, 


C9   =  c2  +  a  C5  "  a 


cio=1-  c3"ac6'   c11=B(c5-i),   cl2=6c6 


Eq.  (46)  can  be  written  in  the  form 


(47) 


K  Au   =  AR 


(48) 
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where 


K  =  K  +  C_  M  +  CQ  KL 


NS 


*St  -  »t+4t  "  M  [C10  ut  -  Cg  ut]  -  K  [„t  -  Cuut  -  C12  Gt]  -  F™      (49) 

Equation  (49)  can  be  solved  for  Au  at  each  time  step  and 
the  nodal  accelerations,  velocities  and  displacements  at  time  t  +  At 
can  be  determined  from  Eqs.  (36)  ~.ud    (44)  giving 


Ht+At  =  Ht  +  Cl  A^t  "  C2  2t  "  S  Ht 


Ht+At  =  2t  +  c4  A^t  "  c5  Ht  "  ce  h  (50) 


u   «   =  u  +  Au 
-t+At   -t    -t 


C.   ACCURACY  OF  SOLUTION 

Equation  (48)  is  an  approximate  form  of  the  actual  equation  of 

motion  to  be  solved  at  each  time  step  since  it  was  derived  from 

approximations  of  the  instantaneous  linearized  stiffness  relations 

between  the  incremental  forces  and  the  incremental  displacements. 

Depending  on  the  nonlinear ity  of  the  equation  and  the  magnitude  of  time 

interval  At,  the  incremental  linearization  may  introduce  an  instability 

in  the  overall  solution.   A  measurement  of  how  well  the  dynamic 

equilibrium  at  time  t+At  is  being  satisfied  by  the  approximate 

solution  of  Eq.  (48)  may  be  expressed  by  the  residual  or  unbalanced 

forces  6r   .  .   The  correctness  of  solution  using  Eq.  (48)  may  then 

be  given  by  comparing  the  ratio  of  the  Euclidean  norm  of  the  residual 

forces  and  the  external  forces  A   with  a  specified  tolerance  A 

P  PS 

using  the  relation 
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where 


&p  =  ll5t+AtIIHl5t+At-^t+Atll  £  ***  <51) 


«?t+At  "  *t+At  "  5  St+it  -  £  St+At  •  5L  2t+At  "  Et+At     <52) 


When  the  accuracy  of  solution  is  unsatisfactory,  it  may  be 
improved  by  using  smaller  time  intervals  or  by  applying  an  equilibrium 
correction  through  an  iteration  process. 

D.   SUBDIVISION  OF  TIME  INTERVAL 

Although  sufficient  accuracy  of  solution  can  be  maintained 
using  small  time  intervals  throughout  the  period  of  dynamic  response, 
the  computational  effort  may  be  excessive.   Therefore,  it  is  more 
efficient  to  use  the  smaller  intervals  only  during  those  intervals  of 
time  when  they  are  actually  needed  and  to  use  longer  intervals  during 
the  remaining  part  of  the  response.   When  the  solution  at  time  t 
satisfies  dynamic  equilibrium  with  the  specified  accuracy  but  the 
solution  at  time  t  +  At  does  not,  one  can  return  to  the  solution  at 
time  t  and  use  smaller  equidistant  sub- intervals  between  times  t  and 
t  +  At.   The  accuracy  of  solution  thus  obtained  by  the  smaller  time 
intervals  can  be  checked  using  Eq.  (51)  at  the  end  of  each  of  the 
smaller  time  intervals. 

As  seen  from  Eq.  (44),  the  nodal  accelerations,  nodal  velocities 
and  nodal  displacements  at  time  t  +  At   can  be  reversed  to  the  solution  at 
time   t  using  the  relations 

Bt  =  C14  A^t  "  C15  Ht+At  "  C16  Ht+At 

2t  =  C17  A^t  "  C18  Ht+At  -  C19  §t+At  (53) 


u  =  u   a   -  Au 
-t    -t+At     -t 
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where 


C  C 

c   =— -i C   -li^  C   --11 

13    6-C7-1/2  '         14    A  2       '  ^15    At 

C   (5-1) 

c16  «  c13(a-6),      c1?  -  — ir—  ,  c18  =  C  (0-1/2)    (54) 


C19  "  C13(f  "  °)At 


E.   EQUILIBRIUM  ITERATION 

The  equation  of  motion  for  the  ith  equilibrium  iteration  at 
time  t  is  expressed  as 


M  6u(1)  +  C  6u(i)  +  K(i)  6u(i)  =  6R(i)  (55) 

-   -t     -   -t     -t    -t      -t 


where,   M  and  C  are  the  constant  mass  and  damping  matrices, 

respectively,   K    is  the  tangent  stiffness  matrix  at  time  t   for  the 

ith  iteration,   6u    ,  6u    ,   and  6u     are  the  corrective  nodal 
-t     -t  -t 

accelerations,  nodal  velocities  and  nodal  displacements,  respectively, 
for  the  ith  iteration  as  defined  by 


...(i)  ...(i+1)  »(i) 

6u  =  u  -  u 

-t  -t  -t 

-•'(i)  '(i+D  *(i)                       ,,,, 

ou  =  u  -  u                             (56) 

*    (i)  (i+D  (i) 

ou  =  u  -  u 

-t  -t  -t 


and  where  6r    representing  the  nodal  residual  forces  for  the  ith 
iteration  is  given  by 


OR    =R-Mu     -Cu     -Ku     _F^  (57) 

-t      -t    -  -t      -  -t      -   -t      -t 
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From  the  integration  scheme  of  Eq.  (44) ,  the  corrective 
accelerations  and  velocities  of  Eq.  (56)  can  be  written  as  [22] 


6u^}  =  c  6u(i) 
-t      1  -t 


(58) 


ou    =  C.  ou 
-t      4  -t 


Substituting  Eqs.  (57)  and  (58)  into  Eq.  (55) ,  one  obtains 


|U>  s^    =  6R^i}  (59) 


where 


K<1}  -  K™    +  C,  M  +  C  pKL 

— t  ~t       /  -     o  — 

OR  =  R  -  M  [u    +  au    ]  -  K  [u     +  pu    ]  -  F           (60) 

-t  -t    -   -t      -t      -   -t      -t      -t 


(i) 


Equation  (59)  can  be  solved  for  the  ith  corrective  displacements  6u 
and  then  the  corrected  nodal  accelerations,  nodal  velocities  and  nodal 
displacements  can  be  calculated  using  Eqs.  (56)  and  (58),  i.e., 

..(i+1)   ..(i)  x  _  .  (i) 

U        =  U      +  C„  OU 

-t     -t     1  -t 

•(i+1)      «(i)  .  (i)  .... 

u      =  u    +  C.  Ou  (61) 

-t       -t      4  -t 

(i+1)     (i)  ^  .  (i) 
u      =  u    +  Ou 
-t       -t      -t 


If  convergence  occurs,  the  iteration  can  be  continued  until 

the  dynamic  equilibrium  of  the  motion  is  satisfied  within  the  specified 

accuracy.   The  convergence  of  equilibrium  iteration  may  be  expressed  by 

comparing  the  ratio  €f  the  Euclidean  norm  of  the  corrective  displacement 

to  the  total  displacement  A   with  a  specified  tolerance  A    in  the  form 

u  us 
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llSt  II  +  llHt    I 


or  by  comparing  the  ratio  of  the  Euclidean  norm  of  the  residual  forces 

to  the  total  external  forces  A   with  a  specified  tolerance  A    in 

p  ps 

the  form 

t  11^"  1 1  x 

A1  =   ^ — <  A1  (63) 

p   Mr  ||  +  ||r  -  6r(i)||  =   ps 

I  l_tM     I  l_t      _t    I  I 

The  step-by-step  integration  algorithm  used  herein  for  non- 
linear systems  including  the  equilibrium  iteration  and  the  subdivision 
of  time  intervals  may  be  summarized  as  follows : 

(1)  Initial  Calculations 

(a)  Form  the  initial  stiffness  matrix  K   and  the  mass  matrix  M 
for  the  system. 

(b)  Solve  for  the  initial  displacement  and  calculate  the  element 
forces  due  to  static  loads. 

(c)  Set  up  the  dynamic  load  and  ground  excitation  time  history. 

•  NS 

(d)  Set   initial   conditions      u     =   u     =  u     =   0,    f      (j)    =0      for 

j£NEL,      Kq  =   KL. 

(e)  Calculate  the  step-by-step  integration  constants  C. 
(i  =  1,2,.... 12)  and  subdivision  control  constants 
C.  (i  =  13,14, 19) . 

(2)  For  Each  Time  Increment  At 
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(f)  Form  the  effective  dynamic  stiffness  matrix  K   and  the 
effective  dynamic  load  vector  Ar  . 

(g)  Solve  for  the  displacement  increment  Au  . 

(h)   Compute  current  accelerations,  velocities  and  displacements 

-t+At'  -t+At'  -t+At" 

(i)   Calculate  current  element  forces,  check  nonlinearity  conditions, 

N 
and  compute  the  new  element  tangent  stiffness  matrix  k  ...   for 

-t(;j) 

j£NEL;  determine  the  inelastic  deformation  vectors,  and  current 

inelastic  element  nodal  restoring  force  vectors   f  ...   for  j^NEL. 

-t(j) 

(j)   Compute  residual  force  vector  6r   .    and  check  the  accuracy 

of  solution; 

i)   if  calculated  responses  are  sufficiently  accurate,  skip 
to  (k) . 

ii)   if  subdivision  of  time  interval  is  needed,  go  to  (3) 

after  reversing  the  current  inelastic  element  deformations 
from  time  t  +  At  to  time  t. 

iii)   if  equilibrium  iteration  is  needed,  go  to  (4) . 

(k)   Repeat  steps  (f)  through  (j)  for  the  next  time  increment. 

(3)   Subdivision  of  Time  Interval 

(£)   Backspace  accelerations,  velocities  and  displacements  from 
time  t+At  to  time  t,   and  calculate  element  tangent 
stiffness  matrices  and  inelastic  element  nodal  restoring 
forces  at  time  t. 

(m)   Calculate  the  dynamic  load  at  time  t  =  t  +  i  At  for 

i=l,2,....n  in  which  At  =  At/n  and  n  is  a  specified 
number  of  time  interval. 
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(n)   Form  the  effective  dynamic  stiffness  matrix  K>  .~  and  the 
effective  dynamic  load  vector  Ar~   ~  . 

(o)   Solve  for  the  displacement  increment  Au~  .~  . 

-t-At 
(p)   Compute  current  accelerations ,  velocities  and  displacements 
u^ ,      u~ ,   u^  . 

-t'  -t'  -t 

(q)   Calculate  current  element  forces ,  check  nonlinearity  conditions , 

N 
and  compute  new  element  tangent  stiffness  matrices  k^  , . .   for 

-t(]) 

j€NEL;  determine  the  inelastic  deformation  vectors,  and  current 

SN 
inelastic  element  nodal  restoring  force  vector  f~  . . .   for 

-t(]) 

jGNEL. 

(r)   Compute  residual  force  vector  6Rp  and  check  the  accuracy  of 

solution; 

i)   if  calculated  responses  are  sufficiently  accurate,  skip 
to  (s)  . 

ii)   if  equilibrium  iteration  is  needed,  go  to  (4) . 

(s)   Repeat  steps  (m)  through  (r)  for  next  sub-time  increment 
(i  =  1,2, . . . .n) .   If  i=n,   then  go  back  to  (k) 

(4)   Equilibrium  Iteration 

(t)   Form  the  effective  dynamic  stiffness  matrix  K   .    and  solve 
for  the  corrective  displacement  vector  6u^_  .   . 

(u)   Compute  the  i+1   corrected  accelerations,  velocities  and 

-.   ,      _   "(i+U    *(i+D     (i+D 
displacements  ufc+At  ,   u^  ,   ufc+At  . 

(v)   Calculate  current  element  forces,  check  nonlinearity  con- 
ditions, and  compute  element  tangent  stiffness  matrices; 
determine  the  inelastic  deformation  vectors,  and  current 
element  nodal  restoring  force  vectors. 
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(w)   Compute  residual  force  vector  6r  .    and  check  the  accuracy 
of  solution; 

i)   if  corrected  responses  are  sufficiently  accurate,  go  back 
to  (k). 

ii)   if  convergence  occurs  and  corrected  responses  are  not  yet 
sufficiently  accurate,  repeat  steps  (t)  through  (w) . 
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IV.   EXPERIMENTAL  MODEL  STUDY  AND  TEST  RESULTS 

The  purpose  of  this  chapter  is  to  present  the  characteristics  of  the 
model  bridge  structure,  test  procedures  and  test  results  obtained  in  the 
Series  IV  excitation  test  program  including  the  results  of  a  number  of 
preliminary  component  tests  used  in  formulating  the  analytical  model  and 
in  interpreting  the  correlation  results  of  seismic  responses  with 
analytical  predictions .   Detailed  information  on  the  experimental  phase 
of  the  overall  program  are  given  in  references  [23]  and  [24] . 

A.   EXPERIMENTAL  MODEL  [23,24] 

The  bridge  model  was  originally  based  on  a  symmetrical  simplified 
version  of  the  east  half  of  the  5/14  South  Connector  Overcrossing  that 
suffered  severe  structural  damages  during  the  San  Fernando  earthquake  in 
1971.   It  was  designed  with  the  aim  of  providing  the  primary  features  of 
high  curved  highway  bridge  structures,  namely,  deck  curvature,  expansion 
joints,  and  long  columns. 

Based  on  size  and  performance  of  the  shaking  table  at  the  University 

of  California,  Berkeley  [27],  scales  of  the  model  were  established  as 

follows; 

L   =   30 
r 

F   =  L2  =  900  (64) 

r      r 

T   =  Sl~     =   1/5.5 
r       r 

in  which  L  ,  F   and  T   represent  the  geometric,  force  and  time  scale 

factors.   In  order  to  bring  the  inertia  forces  including  gravitational 

forces  into  the  force  scale  factor  F  ,   substantial  amount  of  lead 

r 

weights  were  placed  on  the  model.   The  total  weight  including  the  lead 
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weights  was  approximately  8.5  kips  which  corresponds  to  approximately 
8.5%  of  the  weight  of  the  shaking  table. 

The  bridge  model  structure  adopted  was  composed  of  three  sub- 
assemblages  as  shown  in  Fig.  4.1;  a  center  girder/column  system  and  two 
side  girder/column/abutment  systems  which  were  essentially  independent 
of  each  other  having  different  dynamic  characteristics.   These  sub- 
assemblages  were  tied  together  by  two  expansion  joints  placed  in 
symmetrical  positions. 

Seismic  responses  of  displacement  relative  to  the  shaking  table 
were  measured  by  LVDT  at  three  points  of  the  model  as  shown  in  Fig.  4.2, 
i.e.  at  the  end  of  side  girder  No.  1,  center  of  the  center  girder  and 
end  of  side  girder  No.  2.   Two  horizontal  components,  namely  the 
longitudinal  (X)  component  and  the  transverse  (Y)  component,  were 
measured  at  each  point.   In  addition  to  these  displacements,  relative 
movements  between  the  two  end  diaphragms  (opening  and  closing  of  joint 
gap)  of  each  expansion  joint  at  both  the  inner  and  outer  sides  were 
measured.   Thus  in  all,  10  responses  of  displacement  (6  relative  to  the 
shaking  table  and  4  across  the  expansion  joints)  were  recorded  during 
the  periods  of  seismic  excitation.   However  since  the  response  at  the 
center  of  the  center  girder  in  the  longitudinal  (X)  direction  was  very 
small  due  to  the  direction  of  excitation,  it  was  finally  disregarded, 
leaving  9  responses  of  displacement  to  be  used  for  the  correlation  study. 

The  model  was  made  as  a  component  system  allowing  repeated  experi- 
mental data  to  be  derived  from  one  basic  model  as  damaged  components 
were  either  repaired  or  replaced.   The  properties  of  each  component  are 
discribed  as  follows : 

1.   Superstructure  -  The  superstructure  of  the  model  consisted  of 
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one  center  girder  and  two  side  girders.   They  were  constructed  with  a 
rectangular  cross  section  8.5  inches  in  width  and  2  1/2  inches  in  thick- 
ness.  This  section  was  overdesigned  by  a  factor  of  approximately  2  to 
insure  that  the  superstructure  would  behave  elastically  during  all  tests 
[23].   A  high  strength  (8000  psi.)  shrinkage  resistant,  microconcrete  was 
used  in  constructing  the  girders  to  avoid  shrinkage  cracks  which  would 
adversely  affect  the  stiffness  and  damping  characteristics. 

2.  Columns  -  Two  types  of  columns  having  different  sections  were 
used  in  the  model;  i.e.  two  weak  columns  (COL-IB  and  2B)  having  a  section 
2  1/2  x  1  1/2  inches  which  supported  the  center  girder  and  two  strong 
columns  (C0L-1A  and  2A)  having  a  section  2  1/2  x  4  inches  which  supported 
the  side  girders.   All  columns  were  constructed  using  4500  psi.  strength 
normal  Portland  cement  microconcrete  and  four  #2  (2/8  inch)  annealed 
reinforcing  bars  having  a  yield  strength  of  approximately  50  ksi. 

The  bottoms  of  the  columns  were  clamped  to  a  base  steel  frame 
which  was  designed  with  sufficiently  high  stiffness  that  essentially 
full  fixty  was  provided  at  the  base  of  each  column.   The  top  of  the 
columns  were  attached  to  the  girders. 

3.  Expansion  Joints  -  Expansion  joint  element  consisted  of  a  hinge 
seat,  a  transverse  shear  key,  a  rubber  pad,  a  vertical  restrainer  and  a 
pair  of  longitudinal  joint  restrainer  tie  bars  as  shown  in  Fig.  4.3. 

The  hinge  seat  and  the  shear  key  were  constructed  of  concrete 
having  the  same  properties  as  those  of  the  concrete  used  for  the  super- 
structure and  they  were  heavily  reinforced  to  provide  high  resistance  tc 
shear  forces  and  bending  moments.   The  section  of  the  hinge  seat  was 
8.5  inches  in  width  and  0.5  inches  in  length.   The  shear  key  having  a 
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section  2.5  inches  in  width,  0.5  inches  in  length,  and  1.5  inches  in 
height  was  designed  to  prevent  transverse  relative  motion  between  the 
girders. 

A  1/16  inch  thickness  rubber  pad  was  placed  over  the  entire  hinge 
seat  simulating  the  prototype  elastomeric  bearing  pad. 

A  pair  of  vertical  stoppers  was  provided  at  each  expansion  joint  to 
prevent  up-lift  of  the  girders  from  the  hinge  seats. 

A  pair  of  longitudinal  joint  restrainer  bars  5  1/2  inches  in  length 
and  1/8  of  an  inch  in  diameter  was  used  to  tie  together  the  superstructure. 
They  were  made  of  annealed  mild  steels  having  yield  and  ultimate  strengths 
(strain)  of  approximately  650  lbs.  (0.2%)  and  750  lbs.  (20%) ,  respec- 
tively.  The  restrainer  tie  bars  were  mounted  on  each  side  of  the  super- 
structure parallel  to  the  bridge  axis.   They  were  designed  to  resist  only 
opening  of  the  joint  gaps  when  the  relative  movement  exceeded  the 
initial  gap  A  . 

4.   Abutments  -  The  abutments  were  constructed  as  steel  frames  to 
support  the  outer  ends  of  the  side  girders.   They  were  originally 
designed  to  provide  full  fixity  to  the  ends  of  the  girders  except  for 
allowing  rotational  freedom  around  the  transverse  horizontal  (y)  axis. 
However,  due  to  lack  of  stiffness  of  the  steel  cover  plate  at  the  abut- 
ment, some  rotational  freedom  around  the  vertical  (z)  axis  was  also 
provided. 

B.   TEST  PROCEDURES  AND  TEST  RESULTS 

1.   Stiffness  Properties  of  Components  -  A  number  of  preliminary 
tests  were  performed  on  the  components  of  the  model  to  determine  their 
stiffness  properties  and  their  boundary  conditions.   The  test  results 
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obtained  are  used  subsequently  in  formulating  the  analytical  model  of 
the  complete  structure. 

Columns  -  Flexural  rigidities  of  columns  were  determined  from  the 
results  of  static  load-deflection  tests.   The  measured  load  vs.  deflection 
relations  were  linear  up  to  those  levels  later  produced  during  th6  high 
intensity  seismic  excitations.   The  measured  flexural  rigidities  were 
confirmed  by  comparing  the  resulting  calculated  lowest  natural  fre- 
quencies of  the  free  standing  columns  with  their  experimental  values, 

Other  column  rigidities  such  as  those  in  the  longitudinal  and  torsional 

i 

directions  were  calculated  using  the  known  sectional  and  material  pro- 
perties. 

Superstructure  -  The  flexural  rigidities  of  the  girders  around  their 
strong  axes  were  determined  from  static  load  vs.  deflection  tests.   The 
load  vs.  deflection  relations  were  linear  up  to  those  levels  developed 
during  high  intensity  seismic  excitations.   The  flexural  rigidities  of 
the  girders  around  their  weak  axes,  however,  could  not  be  determined  in 
the  same  way  since  twisting  associated  with  the  deck  curvature  was 
coupled  with  bending.   Therefore,  these  rigidities  were  determined  from 
free  vibration  tests  of  the  side  subassemblage  as  described  in  the  sub- 
sequent paragraph.   The  longitudinal  and  torsional  rigidities  were  cal- 
culated from  their  sectional  and  material  properties. 

Expansion  Joint  -  The  expansion  joint  is  the  most  complicated 
component  in  the  entire  bridge  model.   This  complexity  results  from  the 
combined  actions  of  the  rubber  pad,  shear  key  and  tie  bars  associated 
with  slippages  and  collisions  which  take  place  between  girders. 

Figure  4.4  shows  the  force  vs.  relative  displacement  relation  of 
the  expansion  joint  in  which  plus  relative  displacement  corresponds  to 
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an  opening  of  the  joint  gap  and  minus  relative  displacement  corresponds 
closing  of  the  joint  gap.   When  the  expansion  joint  undergoes  an  in- 
creasing relative  movement  in  the  positive  direction,  the  rubber  pad 
first  deforms  under  shear  action  providing  resistance  to  the  motion. 
Then,  slippage  takes  place  when  the  applied  force  reaches  the  maximum 
friction  force  which  can  be  developed  on  the  contact  plane  of  the 
expansion  joint  (Fig.  4.4.1).  It  should  be  noted  here  that  the  difference 
between  dynamic  and  static  frictions  is  disregarded  in  this  force  vs. 
displacement  relation.   When  the  positive  relative  movement  reaches  the 
tie  gap  A  ,   the  pair  of  tie  bars  begin  resisting  further  opening  of 
the  joint  gap.   This  resistance  builds  up  linearly  with  joint  separation 
until  the  yield  strength  of  the  tie  bars  is  reached;  then,  yielding 
under  constant  force  takes  place  (Fig.  4.4.2).  When  the  expansion  joint 
undergoes  a  relative  movement  in  the  negative  direction,  the  rubber  pad 
deforms  and  resists  the  motion  in  the  same  manner  described  above  for 
the  positive  direction.   However  the  tie  bars  do  not  resist  motion  in 
this  direction.   If  the  expansion  joint  undergoes  further  negative 

relative  motion  reaching  the  value  of  the  seat  gap  A  ,   collision 

G 

takes  place  between  the  girders. 

Although  the  expansion  joints  have  such  complex  characteristics,  a 
direct  measurement  of  their  constraint  at  the  expansion  joint  was  not 
performed  because  of  difficulties  involved  in  such  an  experimental  test. 
Therefore  the  experimental  data  obtained  to  formulate  the  analytical 
model  of  the  expansion  joint  were  certain  element  properties. 

The  stiffness  and  strength  of  the  tie  bar  were  measured  after  they 
were  made  by  annealing  mild  steel  [23] .   One  typical  stress  vs.  strain 
relation  is  displayed  in  Fig.  4.5.   A  longitudinal  stiffness  of 
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4 
6.47  x  10  lbs/inch  and  a  yield  (ultimate)  strength  of  650  lbs  (750  lbs) 

were  determined  from  the  test  results. 

Preliminary  tests  were  performed  on  the  rubber  pad  to  determine  its 

shearing  stiffness  [25].   To  accomplish  this,  a  specimen  of  the  rubber 

pad  was  sandwitched  between  two  steel  plates  which  were  then  displaced 

so  that  pure  shear  deformation  was  produced  in  the  pad.   The  specimen 

used  for  this  test  was  1  inch  long  and  1/2  inch  wide.   The  test  was 

performed  with  a  changing  frequency  of  cyclic  shear  load  and  contact 

force.   Some  typical  results  of  the  shear  force  vs.  deformation  are 

shown  in  Fig.  4.6.   From  the  results,  it  is  apparent  that  the  rubber 

pad  has  a  slight  hysteretic  type  behavior  with  the  average  stiffness 

being  significantly  affected  by  contact  pressure  but  less  affected  by 

frequency  in  the  range  0-5  Hz.   Assuming  the  stress-strain  relation  of 

the  rubber  pad  can  be  modeled  by  a  linear  elastic  spring  with  stiffness 

proportional  to  both  the  contact  force  and  the  size  of  specimen,  the 

p 

effective  stiffness  of  the  pad  k   was  found  to  be 

x 

kR  =  22800  lbs/inch  (65) 

x 

A  preliminary  tests  was  also  performed  to  determine  the  friction 
coefficient  [25] .   One  of  the  expansion  joint  components  was  subjected 
to  an  increasing  static  force  in  its  longitudinal  direction  while  under 
compressive  contact  pressure.   The  maximum  frictional  resisting  force 
was  measured  when  slippage  first  occured.   Then,  the  friction  coefficient 
was  obtained  by  dividing  the  measured  maximum  resisting  force  by  the 
contact  pressure.   However,  due  to  large  variations  in  the  measured 
data,  the  test  results  could  only  suggest  that  the  friction  coefficient 
was  within  the  range 
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EJ 
0.3  <  V   <  0.6  (66) 


Therefore  considering  the  condition  of  the  contact  plane,  the  friction 
coefficient  was  taken  as  0.4  for  the  analytical  study.   The  effect  of 
this  assumed  value  on  dynamic  response  is  discussed  in  Chapter  VII. 

Finally,  combining  the  measured  and  estimated  data  for  the  expansion 
joint,  its  force  vs.  relative  displacement  relation  was  established  as 
shown  in  Fig.  4.4.   Since  the  vertical  contact  force  at  the  expansion 
joint  in  the  bridge  model  is  estimated  to  be  approximately  380  lbs,  a 
friction  force  of  approximately  0.4  x  380  lbs  in  the  longitudinal 
direction  can  be  developed.   Under  this  force  the  rubber  pad  can  deform 
approximately  0.01  inch  (0.4  x  380/22,800)  before  slippages  take  place. 
The  maximum  elastic  elongation  of  the  tie  bar  is  approximately  0.01  inch 
(650/64,700)  and  the  maximum  total  tie  bar  force  (2  bars)  is  approximatly 
2  x  650  lbs. 

2.   Static  and  Dynamic  Tests  for  Side  Subassemblage  -  After  one  side 
subassemblage  composed  of  the  side  girder,  column  and  abutment  was 
fablicated,  both  a  free  vibration  test  and  a  static  load-deflection  test 
were  conducted  [25] .   These  test  results  provided  valuable  information  in 
determining  stiffness  properties  of  the  components  and  also  in  formulating 
an  analytical  model. 

The  side  subassemblage  was  first  excited  manually  in  both  horizontal 
and  vertical  directions.   Since  the  frequencies  decreased  slightly  with 
amplitude  of  vibration,  the  frequencies  corresponding  to  those  amplitudes 
developed  during  high  intensity  seismic  excitations  were  estimated 
through  extrapolation  to  give 
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s 

F       -     5.8  Hz       horizontal  vibration 
H 


s 

F   -  9.5  Hz   vertical  vibration 


(67) 


The  side  subassemblage  was  then  subjected  to  an  increasing  static  load 
up  to  400  lbs  in  the  horizontal  radial  direction  and  the  resulting  load 
vs.  deflection  relation  was  measured. 

From  the  test  results,  the  flexural  rigidity  of  the  side  girder 
around  its  weak  axis  was  calculated  from  the  measured  natural  frequency 
of  vertical  vibration.   Since  the  boundary  condition  of  the  side  girder 
in  that  direction  was  well  defined,  the  flexural  rigidity  thus  determined 
was  considered  to  be  sufficiently  reliable. 

The  natural  frequencies  of  the  side  subassemblage  were  then  cal- 
culated based  on  the  known  stiffness  properties  and  boundary  conditions 
and  they  were  compared  with  the  measured  results.   From  these  comparisons 
it  was  found  that  the  calculated  natural  frequency  in  the  horizontal 
direction  was  approximately  25%  higher  than  the  measured  result. 
Because  of  this  discrepancy, the  boundary  condition  of  the  experimental 
model  was  re-examined  and  it  was  found  that  significant  flexibility 
existed  in  the  abutment  allowing  a  rotation  to  develop  around  the 
vertical  (z)  axis  [25] .   The  abutment  steel  cover  was  not  sufficiently 
stiff  to  provide  full  fixity  to  the  side  girder  in  that  direction.   Some 
flexibility  for  the  rotational  component  was  then  introduced  into  the 
analytical  model  so  that  the  calculated  natural  frequency  would  match 
the  measured  value. 

Finally,  a  static  load  was  applied  to  the  corrected  analytical 
model  described  above  and  its  static  load  vs.  deflection  relation  was 
compared  with  the  measured  relations.   As  shown  in  Fig.  4.7,  both  rela- 
tions are  in  good  agreement. 
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Based  on  the  good  agreement  finally  achieved  in  the  static  load  vs. 
deflection  relation,  it  was  felt  that  the  stiffness  properties  and 
boundary  conditions  of  the  side  subassemblage  had  been  determined  with 
sufficient  accuracy  to  be  used  in  formulating  the  analytical  model  of 
the  whole  bridge  model. 

3.   Small  Amplitude  Dynamic  Properties  of  the  Model  -  Free  vibration 
tests  and  forced  harmonic  excitation  tests  were  conducted  to  determine 
the  small  amplitude  dynamic  characteristics  of  the  complete  model.   The 
free  vibration  tests  were  initiated  by  hand  and  it  was  found  that  the 
first  vibration  mode  consisted  of  longitudinal  motion  of  the  center 
girder/column  subassemblage  combined  with  antisymmetrical  motion  of  the 
side  subassemblages  [23] .   The  natural  frequency  of  this  mode  was  found 
to  be 

F   -   5  Hz  (68) 

The  second  free  vibration  mode  was  found  to  be  essentially  a  symmetric 
rigid  body  motion  of  the  superstructure  in  a  horizontal  plane.   The 
natural  frequency  of  this  mode  could  not  be  clearily  obtained  however 
since  the  first  free  vibrations  mode  was  always  excited  making  it 
difficult  to  measure  the  second  mode  frequency.   This  difficulty  was 
encountered  even  when  the  test  was  performed  by  displacing  the  structure 
into  essentially  a  second  mode  shape  and  suddenly  releasing  it  by 
fracturing  the  tensioned  cables  attached  at  three  symmetrical  points  of 
the  model  [23].   It  should  be  noted  that  to  excite  a  particular  mode, 
the  excitation  should  be  applied  so  that  a  maximum  of  energy  is  trans- 
ferred to  this  mode  with  a  minimum  of  energy  transferred  to  its 
neighboring  modes  [33].   As  described  above,  the  vibration  shape  of  the 
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side  subassemblage  were  similar  for  both  the  first  and  the  second  modes; 
thus  increasing  the  difficulty  of  manually  exciting  either  mode  without 
the  other  being  presented.   Furthermore,  the  second  natural  frequency  is 
close  to  the  first  natural  frequency;  thus,  causing  additional  difficulty 
in  distinguishing  the  higher  mode  from  the  lower  mode. 

Forced  harmonic  excitation  tests  were  conducted  by  sweeping  the 
frequency  of  table  input  motion  from  5.5  Hz  to  7.5  Hz  for  transverse 
vibration  and  from  8  Hz  to  10  Hz  for  vertical  vibration.   The  table 
acceleration  amplitude  was  held  constant  over  both  frequency  ranges. 
From  these  tests ,  it  was  found  that  the  transverse  vibration  mode  was 
essentially  the  same  as  that  observed  from  the  free  vibration  test  and 
the  vertical  vibration  mode  produced  vertical  motions  primarily  in  the 
side  girders.   Resonant  curves  of  response  amplitude  are  shown  in 
Fig.  4.8  for  the  transverse  vibration.   The  natural  frequency  and 
damping  ratio  determined  from  this  test  were 


Fm  -   6.6  Hz 
T 


£   -   6%  of  critical 


transverse  vibration 


(69) 


For  vertical  vibrations,  the  amplitudes  could  not  be  easily  measured  so 
that  only  the  natural  frequency  was  measured,  i.e., 


F   -  9     10  Hz   vertical  vibration 


(70) 


4.   Seismic  Behavior  of  the  Model  -  The  model  was  subjected  to  six 
different  earthquake  excitations  in  the  Series  IV  test  program.   Three 
of  these  excitations  designated  as  Tests  Hi,  H2 ,  and  H3  consisted  of 
horizontal  motion  only  while  the  other  three  excitations  designated  as 
Tests  HV1,  HV2,  and  HV3  consisted  of  both  horizontal  and  vertical  motions, 
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The  horizontal  excitation  was  prescribed  by  an  artificially  generated 
accelerogram  previously  defined  for  post  earthquake  studies  of  the  Olive 
View  Hospital  [28] .   The  vertical  excitation  was  also  prescribed  by  an 
artificially  generated  accelerogram  [24].   The  peak  vertical  excitation 
was  chosen  to  be  approximately  one  half  the  peak  horizontal  excitation. 
Both  vertical  and  horizontal  excitations  were  scaled  in  time  1/5.5  so  as 
to  satisfy  the  time  scale  factor  T  .   The  acceleration  response  spectra 
for  the  motions  measured  on  the  shaking  table  during  Tests  Hi,  H3,  and 
HV2  are  shown  in  Fig.  4.9.   It  is  observed  that  although  there  are  some 
discrepancies  in  the  amplification  factors,  the  maximum  response 
spectral  values  occured  at  approximately  10  Hz  for  both  the  horizontal 
and  vertical  motions. 

Peak  table  accelerations,  maximum  horizontal  response  of  center 
girder,  initial  tie  gaps,  and  observed  damages  for  the  six  excitation 
tests  are  summarized  in  Table  4.1.   The  superstructure  and  the  abutments 
suffered  no  visiable  structural  damages  throughout  this  series  of  tests. 
The  columns  suffered  no  damage  except  for  one  weak  column  that  suffered 
a  slight  hair  crack  at  its  base.   The  shear  keys  of  the  expansion  joints, 
especially  at  expansion  joint  No.  2,  suffered  severe  spalling  of  concrete 
during  Tests  H3  and  HV3  requiring  that  they  be  replaced  after  these 
particular  tests.   The  tie  gaps  were  set  at  almost  zero  values  before 
Tests  Hi  and  HV1 ;  however,  for  the  other  tests,  the  initial  tie  gaps  were 
controlled  by  the  previously  accumulated  elongations  of  the  tie  bars. 
Collisions  occurred  in  the  expansion  joints  during  Test  H3,  HV1,  HV2  and 
HV3.   Sometimes  during  high  intensity  excitations,  impact  at  one  joint 
would  immediately  cause  impact  to  occur  at  the  other  joint.   Thus, 
multiple  collisions  were  observed  between  the  center  girder  and  the  two 
side  girders.   The  sounds  produced  by  the  impacts  were  quite  severe  in 

intensity. 
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The  relation  between  peak  horizontal  table  acceleration  and  maximum 
horizontal  displacement  at  the  deck  center  is  plotted  in  Fig.  4.10. 
From  this  relation,  it  is  observed  that  the  maximum  response  of  the  model 
increases  almost  linearly  with  peak  table  acceleration  in  the  lower 
intensity  range;  however,  a  scatter  of  maximum  response  values  is 
observed  in  the  higher  intensity  range,  i.e.  for  Tests  H3 ,  HV2  and  HV3. 
Although  Tests  H3  and  HV2  were  conducted  using  almost  the  same  intensity 
of  horizontal  table  accelerations,  namely  0.5G,  the  maximum  response  for 
Test  H3  is  approximately  50%  larger  than  that  produced  in  Test  HV2. 
Certainly  the  presence  of  the  vertical  excitation  in  Test  HV2  cannot 
be  credited  with  this  50%  difference.   The  same  discrepancy  is  observed 
when  comparing  the  results  obtained  in  Tests  HV2  and  HV3.   Test  HV3 
followed  immediately  after  the  Test  HV2  and  used  almost  the  same 
intensities  of  table  motion.   To  explain  this  discrepancy, it  is 
necessary  to  recall  that  a  shear  key  in  expansion  joint  No.  2  which  was 
undamaged  after  Test  HV2  suffered  severe  spalling  failure  during  Tests  H3 
and  HV3.   Thus,  it  is  believed  that  failure  of  the  shear  key  allowed 
large  transverse  relative  motions  to  take  place  between  the  girders  at 
the  expansion  joint,  which,  in  turn, resulted  in  a  significant  increase 
in  overall  response  of  the  complete  bridge  model.   It  is  of  interest  to 
speculate  why  the  shear  key  failed  during  Tests  H3  and  HV3  but  did  not 
in  the  previous  Test  HV2.   The  initial  and  boundary  conditions  for  the 
superstructure,  columns  and  abutments  are  the  same  for  all  three  tests; 
however,  the  initial  tie  gaps  were  different.   For  Test  HV2  the  initial 
tie  gaps  were  small  since  little  yielding  of  the  tie  bars  had  taken  place 
during  Test  HV1.   Due  to  yielding  of  the  tie  bars  which  occurred  during 
the  previous  tests,  the  initial  gaps  were  reasonably  large  for  Tests  H3 
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and  HV3.   It  should  be  noted  that  when  superstructure  motion  takes  place 
in  the  outward  direction,  only  the  tie  bars  resist  the  motions  at  the 
expansion  joints.   Hence  the  large  initial  tie  gaps  in  Tests  H3  and  HV3 
allowed  the  model  to  build  up  large  transverse  oscillations  before  the 
tie  bars  became  effective.   These  large  oscillations  were  sufficient  to 
cause  failure  of  the  shear  key. 

Based  on  the  above  considerations  of  model  response ,  three  tests 
were  chosen  among  the  six  conducted  for  correlation  purposes,  namely, 
Tests  Hi,  H3  and  HV2.   Test  HI  represents  a  typical  response  for  low 
intensity  excitation.   No  yielding  of  tie  bars  or  collisions  of  girders 
took  place  in  this  test.   Tests  H3  and  HV2  on  the  other  hand  represent 
typical  response  for  high  intensity  excitations.   Yielding  of  the  tie 
bars  and  multiple  collision  did  take  place  in  these  tests.   In  addition, 
a  failure  of  the  shear  key  is  believed  to  have  taken  during  Test  H3. 

Test  Hi  -  The  model  was  subjected  to  low  intensity  excitation  with 
a  peak  acceleration  of  0.11G  in  the  transverse  (Y)  horizontal  direction. 
The  table  acceleration  and  the  measured  responses  of  displacement  are 
displayed  in  Figs.  4.11  and  4.12,  respectively. 

During  this  test,  collisions  of  the  girders  did  not  occur  since  the 
relative  response  displacements  at  the  expansion  joints  remained  below 
the  initial  joint  gap  of  approximately  0.05  inch.   The  joint  restrainer 
tie  bars  however  resisted  the  joint  separations.   It  was  found  by 
inspection  following  the  test  that  yielding  of  the  tie  bars  had  not  taken 
place.   Since  the  tie  bars  resisted  only  joint  separation,  the  displace- 
ment response  was  small  in  the  outward  direction  but  large  for  the  inward 
direction.   While  some  non-symmetrical  model  response  was  observed,  the 
major  response  was  primarily  in  a  symmetrical  mode.   This  type  of  response 
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should  be  expected  due  to  the  symmetry  conditions  of  model  and  seismic 
loading. 

Test  H3  -  During  this  test,  the  model  was  subjected  to  high  intensity 
excitation  having  a  peak  acceleration  of  0.5G  in  transverse  (Y)  horizontal 
direction.   The  table  acceleration  and  the  measured  responses  of  dis- 
placement are  displayed  in  Figs  4.13  and  4.14,  respectively. 

Multiple  collisions  of  the  girders  together  with  yieldings  of  the 
joint  restrainer  tie  bars  took  place  at  the  expansion  joints.   Further- 
more ,  spalling  failure  of  the  transverse  shear  key  took  place  at 
expansion  joint  No.  2.   Like  Test  Hi,  the  responses  of  displacement, 
especially  at  the  expansion  joints,  were  somewhat  unsymmetrical.   How- 
ever they  were  of  a  type  opposite  to  that  developed  in  Test  HI,  i.e.  the 
motions  were  large  in  the  outward  direction  and  small  in  the  inward 
direction.   This  type  of  response  was  caused  by  closure  of  the  joints 
during  inward  motion  causing  stiff  arch  action  to  take  place  between 
abutments  while  motion  in  the  outward  direction  was  resisted  mainly  by 
the  more  flexible  tie  bars.   The  maximum  opening  at  expansion  joint 
No.  2  equal  to  approximately  0.35  inch  exceeded  that  of  expansion  joint 
No.  1  by  a  factor  of  approximately  2.   Hence  adding  the  maximum  closure 
(initial  joint  gap)  of  approximately  0.05  inch  to  the  maximum  opening, 
the  total  relative  displacement  of  the  joint  was  approximately  0.4  inch 
which  corresponds  to  4/5  of  the  hinge  seat  length.   This  implies  that 
the  center  girder  was  supported  on  only  a  0.1  inch  hinge  seat  at  the 
expansion  joint.   The  damages  of  the  rubber  pad  observed  after  this  test 
were  undoubtedly  caused  by  this  severe  relative  motion. 

It  was  observed  from  the  measured  response  that  substantial  drift, 
especially  at  expansion  joint  No.  2,  took  place  during  the  excitation. 
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Once  this  drift  developed,  the  subsequent  responses  did  not  return  to 
their  original  values  but  simply  vibrated  the  new  drifted  positions.   It 
should  be  noted  that  although  the  spalling  failure  of  the  transverse 
shear  key  allowed  a  large  relative  transverse  movement  to  take  place 
between  the  girders,  it  had  essentially  nothing  to  do  with  the  observed 
drift  in  response.   Most  likely  this  drift  was  caused  by  a  combination 
of  concrete  fractures  in  the  shear  key  and  tearing  of  the  rubber  pad  in 
the  joint.   Obviously,  such  drift  could  not  be  predicted  analytically. 

Test  HV2  -  During  this  test,  the  model  was  subjected  simultaneously 
to  high  intensity  horizontal  and  vertical  excitations  having  peak 
accelerations  of  0.47G  and  0.27G,  respectively.   The  table  accelerations 
and  the  measured  responses  of  displacement  are  displayed  in  Figs.  4.15 
and  4.16,  respectively.   The  general  features  of  this  test  are  similar 
to  those  of  Test  H3  except  that  failure  of  the  shear  key  did  not  take 
place.   Therefore,  the  overall  response  of  the  model  was  more  symmetrical 
than  the  response  of  Test  H3  and  drifts  in  response  did  not  develop. 
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V.   FORMULATION  OF  A  BASIC  ANALYTICAL  MODEL 

The  purpose  of  this  chapter  is  to  integrate  the  component 
stiffness  properties  and  boundary  conditions  previously  described  into 
a  basic  analytical  model  that  satisfies  the  dynamic  characteristics  of 
the  experimental  model  observed  under  small  amplitude  free  vibration  and 
forced  harmonic  conditions.   These  small  amplitude  characteristics  are 
considered  a  fundamental  requirement  of  the  analytical  model  to  be  used 
in  predicting  seismic  responses  of  the  experimental  model. 

A.  IDEALIZATION  OF  EXPERIMENTAL  MODEL 

During  the  free  vibration  and  forced  harmonic  excitation  tests, 
all  components  of  the  bridge  model  remained  elastic.   Therefore,  the 
bridge  model  was  treated  as  a  linear  three-dimensional  discrete  parameter 
system  consisting  of  55  nodal  points.   In  this  idealization,  24  linear 
curved  beam  elements  respresenting  the  superstructure,  24  linear  straight 
beam  elements  representing  the  columns,  2  linear  elements  representing 
the  expansion  joints  and  2  boundary  elements  representing  the  abutments 
were  used  as  shown  in  Fig.  5.1. 

B.  BASIC  ANALYTICAL  MODEL 

The  component  stiffness  properties  and  boundary  conditions 
previously  presented  are  summarized  again  in  Table  5.1.   Except  for  the 
expansion  joint  properties,  these  characteristics  have  been  fully 
verified  for  including  in  the  basic  analytical  model.   The  various 
constraints  provided  by  the  expansion  joint  need  further  clarification 
however  since  they  are  important  factors  controlling  the  dynamic  response 
characteristics  of  the  analytical  model.   As  previously  shown  in  Fig.  4.4, 
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the  force  vs.  relative  displacement  relation  of  the  expansion  joints  is 
very  complicated  showing  numerous  changes  in  stiffness  with  relative 
displacement.  When  the  relative  displacement  is  very  small,  the  effects 
of  tie  bars  are  negligible  and  slippages  and  collisions  do  not  occur. 
In  such  cases,  only  the  transverse  shear  key,  the  vertical  restrainers 
and  the  rubber  pad  resist  the  relative  movement  between  girders.   There- 
fore in  defining  the  analytical  model,  it  is  necessary  that  these  con- 
straints under  small  amplitude  response  be  adequately  modeled.  First, 
a  set  of  linear  elastic  springs  were  selected  to  resist  the  relative 
translational  and  rotational  movements  between  the  girders.   The 
longitudinal  stiffness  k    and  the  rotational  stiffness  kfi   around 
a  vertical  (z)  axis  were  evaluated  considering  the  shear  stiffness  of 
the  rubber  pad;  thus, 


EJ     R  4 

k   -  k  =  2.28  x  10  lbs/inch 
x     x 


EJ     2        R  5 

kQ   -  d  /12  •  k  =  1.4  x  10  lb.  inches/rad. 
6z  x 


(71) 


in  which  d  represents  width  of  the  superstructure.   Next,  the  vertical 
restrainers  and  the  transverse  shear  key  were  assigned  infinite  stiffness 
to  prevent  vertical  and  transverse  relative  motions,  respectively.  With 
these  prescribed  constraints,  three  degrees  of  freedom  were  permitted  in 
the  joint,  namely,  rotations  around  x  and  y  axes  and  uniform  joint 
separation.   These  prescribed  conditions  (71)  together  with  the  stiffness 
and  boundary  conditions  of  Table  5.1  resulted  in  the  following 
analytically  predicted  frequencies 

f  =  4.1  Hz    longitudinal  vibration 

Li 

fm  =  7.0  Hz    transverse  vibration  (72) 

T 

f  =  9.2  Hz    vertical  vibration 
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As  described  in  Chapter  IV,  the  measured  natural  frequencies  of  the 
experimental  model  were 

F  -  5  Hz         longitudinal  vibration 

L 

F  -  6.6  Hz       transverse  vibration  (73) 

F  -  9  ^  10  Hz    vertical  vibration 

Comparing  both  sets  of  frequencies,  it  is  apparent  that  the  calculated 
frequencies  are  fairly  close  to  their  corresponding  experimental  values . 
Even  so,  this  discrepancy  was  judged  to  be  too  great  to  satisfy  modelling 
of  the  small  amplitude  dynamic  characteristics  of  the  experimental  model. 
The  reason  for  the  discrepancy  derived  from  disregarding  the  effects  of 
tie  bars,  shear  key,  etc* 

To  remove  the  discrepancy,  a  parameter  study  was  carried  out  by 
changing  the  pair  of  expansion  joint  stiffnesses.   First,  this  study 
provided  information  on  the  effects  of  each  individual  stiffness  on  the 
low  frequency  characteristics  of  the  model,  and  second  it  indicated  the 
proper  combination  of  stiffnesses  to  give  best  agreement  of  analytical 

and  experimental  frequencies.   Table  5.2  shows  the  final  results  of  this 

E  J   E  J 
parameter  study.   Clearly,  they  show  that  stiffnesses  k   ,  k_   and 

k    are  more  sensitive  parameters  than  are  stiffnesses  kn   and  kn 
y  9x        6y 

The  effects  of  the  former  three  parameters  on  the  calculated  natural 
frequencies  may  be  written  as 


f   oc  V 

L    0z 


.  E  J     ,   ,  E  J 
f „  x   k    and  k 
T    x        y 


(74) 
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Further,  it  was  found  that  the  combination  of  stiffnesses  giving  best 
correlations  of  frequencies  was 


kEJ  =  1.9  x  104  lbs/inch 

kEJ  =  5  x  104  lbs/inch  (75) 

y 

EJ     7 
kQ  =  10  lb.  inches/rad. 

Uz 


The  resulting  natural  frequencies  were 


f  =  5.0  Hz 


f  =  6.7  Hz  (76) 


f   =  9.3  Hz 


It  is  interesting  to  note  that  the  primary  changes  in  stiffnesses  given 

EJ 

by  Eqs.  (71)  and  (75)  are  the  introduction  of  k    and  a  significant 

EJ  E  T 

change  in  kA   .   Transverse  stiffness  k    which  resists  relative 
6z  y 

transverse  movement  between  the  girders  was  incorporated  into  the 
analytical  model  not  only  to  adjust  the  calculated  natural  frequencies 
but  also  because  observations  during  experimental  tests  indicated  that 
small  transverse  relative  motions  actually  took  place  during  high 
intensity  excitations  eventhough  the  transverse  shear  key  prevented 

large  relative  motions  to  occur  [26] .   The  high  stiffness  of 

4  EJ 

5  x  10   lbs/inch  selected  for  k    was  considered  realistic  in  allowing 

small  transverse  relative  motions  to  develop  in  the  analytical  model. 

The  change  in  stiffness  k„    can  be  explained  since  kn   as  given  by 

DZ  OZ 

Eq.  (71)  was  derived  considering  only  the  shear  stiffness  of  the  rubber 
pad.   Actually,  however,  the  transverse  shear  key  gives  some  contribution 
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E  J 
to  k-    since  its  corners  and  edges  resist  rotational  movement,  and  the 
Wz 

joint  restrainer  tie  bars  also  make  a  contribution;  thus,  the  incremental 

EJ 
change  in  kft    from  Eq.  (71)  to  Eq.  (75)  is  reasonable. 

Finally,  the  calculated  lowest  vibration  mode  shapes  in  the 
three  major  directions,  namely,  transverse,  longitudinal  and  vertical, 
were  checked  against  the  observed  vibration  modes  of  the  experimental 
model.   As  is  shown  in  Fig.  5.2,  it  is  apparent  that  all  three  of  the 
calculated  modes  are  very  similar  to  the  corresponding  experimental 
modes. 

Based  on  the  discussion  above,  it  is  clear  that  the  analytical 
model  having  stiffness  properties  and  boundary  conditions  as  shown  in 
Table  5.1  and  having  the  idealized  expansion  joint  constraints  as 
represented  by  Eq.  (75)  adequately  satisfies  the  small  amplitude  dynamic 
characteristics  of  the  experimental  model.   It  should  be  noted  however 
that  since  the  spring  stiffnesses  representing  the  expansion  joint 
constraints  were  modeled  so  that  the  low  frequency  characteristics  of 
the  analytical  model  correspond  with  those  of  the  experimental  model, 
all  errors  involved  in  determining  the  stiffness  properties  and  the 
boundary  conditions  of  the  model  are  absorbed  in  the  selected  set  of 
spring  stiffnesses  of  expansion  joints. 
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FIG.  5.1   ANALYTICAL  MODEL  OF  BRIDGE  MODEL 


FIG.  5.2.1   CALCULATED  LOWEST  VIBRATION  MODE;  LONGITUDINAL  (X)  DIRECTION 
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FIG.  5.2.2   CALCULATED  LOWEST  VIBRATION  MODE;  TRANSVERSE  (Y)  DIRECTION 


FIG.  5.2.3   CALCULATED  LOWEST  VIBRATION  MODE;  VERTICAL  (Z)  DIRECTION 
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VI.   ANALYTICAL  PREDICTION  OF  SEISMIC  RESPONSES 

The  purpose  of  this  chapter  is  to  correlate  the  measured  responses 
of  displacement  of  the  experimental  model  with  the  linear  and  nonlinear 
responses  determined  analytically  for  the  three  seismic  excitations, 
i.e.  for  Tests  HI,  H3  and  HV2.   The  correlation  was  initiated  with 
Test  HI,  since  this  test  represented  essentially  linear  elastic  behavior. 
The  correlation  derived  from  this  test  was  basic  as  it  also  represented 
the  initial  low  amplitude  behavior  of  the  model  during  Tests  H3,  and 
HV2  which  subsequently  experienced  high  amplitude  nonlinear  behavior. 

The  correlation  was  made  by  comparing  predicted  and  measured 
responses,  paying  particular  attention  to  the  agreement  of  maximum 
responses,  phase  relations,  and  the  overall  dynamic  behavior. 

A.   LOW  INTENSITY  SEISMIC  EXCITATION  (TEST  HI) 

1.   Linear  Analysis  Correlation  -  The  basic  analytical  model  that 
satisfied  the  small  amplitude  dynamic  characteristics  of  the  experimental 
model  was  used  in  predicting  seismic  responses  for  Test  HI.   It  was 
subjected  to  the  table  acceleration  time-history  shown  in  Fig.  4.11  in 
its  transverse  direction.   A  time  interval  of  0.01025  second  was 
adopted  for  the  numerical  integration.  A  damping  ratio  of  6%  of 
critical  was  assumed  based  on  experimental  evidence. Eq.  (69)  obtained 
from  low  intensity  forced  harmonic  tests  at  amplitudes  of  vibration 
comparable  to  those  developed  in  the  Test  HI. 

Displacement  response  time-histories  were  calculated  under  these 
conditions  and  were  compared  with  the  corresponding  experimental 
results.   From  these  comparisons,  it  was  found  that  although  global 
behavior  of  the  predicted  response  was  similar  to  the  measured  response, 
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the  correlation  of  individual  peak  amplitudes  and  their  phase  relations 
was  relatively  poor.   Further,  the  predominant  frequency  of  the  pre- 
dicted response  was  slightly  higher  than  that  of  the  measured  response. 
This  implied  that  the  total  stiffness  of  the  basic  analytical  model  was 
slightly  higher  than  that  of  the  experimental  model  under  this 
particular  seismic  excitation.   It  was  then  decided  to  modify  the  stiff- 
ness of  the  analytical  model  slightly  in  order  to  make  the  phase  lag 
discrepancies  as  small  as  possible.   The  overall  stiffness  of  the 
expansion  joints  was  chosen  for  this  purpose  since  it  dominated  the  low 
frequency  characteristics  of  the  model  and  also  since  it  was  the  most 
uncertain  property  established  from  experimental  evidence.   Based  on 
the  parameter  study  represented  by  Eq.  (74) ,  it  was  apparent  that  either 

E  J        EJ 

stiffness  k    or  k    should  be  changed  in  order  to  reduce  the  natural 
x        y 

frequency  of  vibration  in  the  transverse  direction.   As  described  in 

EJ 
the  previous  chapter,   k    was  introduced  into  the  basic  analytical 

model  to  permit  some  small  transverse  relative  motions  to  develop  at 

EJ 
the  expansion  joints.   It  was  clear  therefore  that  stiffness  of  k 

should  not  be  decreased  from  its  original  value  since  such  a  reduction 

would  cause  large  transverse  relative  motions  to  take  place.   On  the 

EJ 
other  hand,  stiffness  k    was  derived  from  the  combined  effects  of 

x 

the  tie  bars,  slippages,  etc.;  thus,  significant  uncertainty  was  present 

in  setting  its  numerical  value.   Hence,   k    was  chosen  to  adjust  the 

x 

frequency  and  after  some  calculations ,  it  was  set  at  the  value 


P  7 

k   =  18400  lbs/inch  (77) 

x 


giving  lowest  natural  frequencies  as  follows : 
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f  =  4.9  Hz  longitudinal  vibration 

Xj 

f  =  6.6  Hz  transverse  vibration  (78) 

f  =  9.3  Hz  vertical  vibration 


The  resulting  predicted  response  of  displacement  is  shown  in  Fig.  6.1 
together  with  the  measured  result.   From  this  comparison  it  is  seen  that 
the  predicted  response  agrees  fairly  well  with  that  of  the  experimental 
model.   The  amplitudes  of  the  predicted  response  are  very  close  to 
those  of  the  response  which  indicates  that  the  established  damping  ratio 
of  6%  was  appropriate.   As  one  would  expect,  the  predicted  response  is 
essentially  symmetric  both  in  the  outward  and  inward  directions  while 
the  measured  response  is  somewhat  unsymmetric  due  to  differences  in 
the  joint  restrainer  tie  bar  properties.   This  feature  can  be  observed 
more  clearly  by  comparing  the  responses  of  expansion  joints. 

The  above  correlation  clearily  indicates  that  the  measured  response 
during  Test  HI  could  be  predicted  with  fairly  good  accuracy  using  a 
linear  analytical  model  provided  its  fundamental  natural  frequency  in 
the  direction  of  excitation  agreed  well  with  the  experimental  value  and 
provided  the  proper  damping  factor  was  introduced  into  the  analysis. 
It  is  quite  evident  however  that  the  unsymmetric  response  caused  by 
differences  in  the  discontinuous  expansion  joint  constraints  cannot  be 
predicted  using  the  linear  analytical  model. 

2.   Nonlinear  Analysis  Correlation  -  For  this  correlation  study, 

the  linear  expansion  joint  element  previously  used  was  replaced  by  the 

nonlinear  element.   It  was  necessary  therefore  to  change  the  stiffnesses 

associated  with  the  tie  bars.   In  the  linear  analysis  correlation,  the 

constraints  of  the  expansion  joints  were  modeled  by  the  following  set 

of  linear  spring  stiffnesses: 
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FT  F*  T 

k   =  18400 ^lbs/inch,  k   =  50000  lbs/inch 

EJ     7  <79> 

kQ  =10  lb.  mches/rad. 
t)z 


For  the  present  nonlinear  analysis  correlation,  it  is  not  necessary  to 

F  -r 

change  k     since  the  tie  bars  do  not  affect  relative  transverse 

y 

EJ 
motion  at  the  expansion  joints.   On  the  other  hand,  stiffnesses   k 

EJ 

and  k„   are  directly  dependent  upon  the  tie  bars.   Therefore,  they 

oz 

must  be  modified.   To  accomplish  these  changes,  the  force  vs.  relative 
displacement  relation  shown  in  Fig.  4.4  was  used.   In  Test  HI  the 
maximum  relative  movements  developed  at  the  expansion  joints  were 
approximately  +  0.01  inch  (seperation)  and  -  0.025  inch  (closing). 
When  a  Coulomb  friction  coefficient  of  0.4  is  prescribed,  the  elastic 
shear  deformation  of  the  rubber  pad  is  approximately  0.01  inch  beyond 
which  slippage  takes  place  between  the  contact  surface  of  the  expansion 
joint.   By  comparing  the  amount  of  relative  movement  developed  during 
excitation  and  the  amount  of  shear  deformation  possible  in  the  rubber 
pad,  it  is  clear  that  maximum  slippages  of  the  order  of  0.015  inch  took 
place  at  the  expansion  joints.   Since  the  amount  of  slippage  in  Test  Hi 
was  very  small ,  it  was  decided  not  to  introduce  Coulomb  type  friction 
force  in  this  analysis.   It  was  then  decided  to  take  account  of  the 
small  slippages  approximately  by  linearlizing  the  force  vs.  relative 

displacement  relation.   Thus,  the  longitudinal  linearized  stiffness 

EJ 
k   0  was  estimated  from  Fig.  4.4  as 

X  J  A/ 

k   0   -   9000  lbs/inch  (80) 

X  ,  A/ 

EJ 
The  stiffness  k„   was  also  changed  for  the  same  reason  described  above 

DZ 

resulting  in  the  value 
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EJ  6 

k„   =  6  x  10   lb.  inches/rad.  (81) 

uz 


Damping  coefficients  a  and  $  were  assigned  values  so  that  the 
damping  ratio  of  6%  of  critical  would  be  achieved  for  the  first  and 
second  normal  modes  of  vibration  as  calculated  by  the  linear  analytical 
model.   The  initial  tie  gaps  A   were  assigned  zero  values  based  on 
measurements  taken  just  prior  to  the  test  which  showed  that  they  were 
very  small. 

The  nonlinear  analytical  model  defined  above  was  subjected  to  the 
table  acceleration  time-history  as  used  in  the  linear  analysis 
correlation.   A  time  interval  of  0.01025  second  was  adopted  for  the 
numerical  integration.   Equilibrium  iteration  was  used  both  alone  and 
with  the  subdivision  of  time  intervals  into  3  equidistant  sub-time 
intervals,  when  needed.   The  tolerances  used  to  control  the  equilibrium 
iteration  and  the  subdivision  defined  in  Eqs.  (51)  and  (62)  were 


A    =  0.01   equilibrium  iteration 
ps 

A    =0.08   subdivision  (82) 

ps 

A1   =  0.001 
us 


The  predicted  displacement  response  was  compared  with  the  experimentally 
measured  result  and  it  was  found  that  the  predicted  response  was  sub- 
stantially improved  as  compared  with  the  linear  analysis  correlation 
although  discrepancies  were  still  observed  both  in  the  peak  amplitudes 
and  their  phase  relations.   The  largest  discrepancy  was  seen  at  the 
expansion  joints  where  openings  of  the  gaps  were  excessively  sup- 
pressed by  the  tie  bars.   Because  of  this  discrepancy,  it  was  decided 
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to  change  the  originally  assumed  zero  tie  gaps  to  small  realistic 

EJ 

values.  Hence,  both  the  assumed  linearized  stiffness  of  k  „  and 

x,£ 

the  initial  tie  gaps  were  adjusted;  finally,  giving  good  correlation 
under  the  conditions 


EJ 
k  0     =     10000  lbs/inch 

(83) 

A     =  0.05  inch  for  both  expansion  joints 


The  comparative  plots  of  predicted  and  measured  responses  are  shown 
in  Fig.  6.2.   Excellent  agreement  can  be  seen  for  these  responses  at  the 
center  point  of  the  model  and  at  both  expansion  joints,  and  the  un- 
symmetrical  vibration  shown  by  the  experimental  response  is  well  followed 
in  the  analysis.   The  predicted  responses  of  force  induced  in  the  tie 
bars  are  shown  in  Fig.  6.3.   Since  the  maximum  value  of  the  predicted  tie 
bar  force  is  approximately  400  lbs,  yielding  of  the  tie  bar  was  not 
indicated  by  analysis  which  is  consistent  with  the  observed  results. 

B.   HIGH  INTENSITY  SEISMIC  EXCITATION  (TEST  HV2) 

1.   Linear  Analysis  Correlation  -  To  carry  out  the  linear  analysis 
correlation  for  Test  HV2,  the  same  linear  model  previously  described  for 
Test  HI  was  used.   This  analytical  model  was  subjected  to  the  measured 
table  accelerations  in  both  the  transverse  and  vertical  directions.   A 
time  interval  of  0.01026  second  was  used  throughout  the  correlation. 
Although  it  was  somewhat  questionable  if  the  same  viscous  damping  ratio 
as  used  in  Test  HI  should  be  used  for  this  test  due  to  the  high  amplitude 
motion,  it  was  decided  to  use  the  same  6%  of  critical  as  a  first  trial. 

The  responses  predicted  on  this  basis  were  compared  with  the 
measured  results  as  shown  in  Fig.  6.4.   As  can  be  seen,  the  predicted 
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response  is  significantly  different  from  the  measured  response  both  in 
the  magnitude  of  peak  amplitudes  and  in  the  frequency  characteristics . 
It  soon  became  apparent  from  this  correlation  that  several  collisions 
within  the  joints  and  yielding  of  the  tie  bars  do  not  allow  the  linear 
analysis  acceptable  correlations.   Therefore,  to  obtain  reasonable 
correlation,  one  must  introduce  these  effects  into  the  analytical  model. 
Thus,  further  correlation  using  the  lienar  analytical  model  after 
adjusting  parameters  was  not  attempted  for  Test  HV2. 

2.   Nonlinear  Analysis  Correlation  -  In  order  to  introduce  the 
effects  of  collisions  and  yielding  of  tie  bars  into  the  analysis,  the 
nonlinear  analytical  model  formulated  for  the  correlation  of  the  Test 
Hi  was  further  modified. 

First,  the  modeling  of  constraints  in  the  expansion  joint  was 
reviewed.   The  assumed  force  vs. relative  displacement  relation  shown  in 
Fig.  4.4  was  again  used  for  this  purpose.   A  principal  difference  in  the 
seismic  responses  between  this  test  and  the  Test  HI  is  the  amplitude 
level.   For  example,  the  maximum  relative  displacements  at  expansion 
joint  No.  2  are  approximately  +  0.2  inch  (opening)  and  -  0.05  inch 
(closing)  in  this  test  which  are  roughly  10  times  those  values  developed 
during  Test  Hi.   It  is  apparent  from  these  large  relative  movements 

that  slippages  of  large  amplitude  took  place  during  the  excitation. 

EJ 
Thus,  instead  of  using  the  linearized  stiffness  k  0     as  in  Test  HI, 

one  must  now  introduce  the  Coulomb  friction  which  can  be  modeled  by  an 

elastoplastic  hysteretic  relation.   The  initial  elastic  slope  of  this 

Coulomb  friction  relation  was  taken  equal  to  the  shear  stiffness  of  the 

rubber  pad  k   .   By  doing  so,  the  assumed  relation  of  Fig.  4.4  can  be 

ri-r  PT 

exactly  modeled  in  the  analysis.   Stiffnesses  k    and  k„   were 
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unchanged  from  Test  HI  since  they  do  not  affect  the  Coulomb  friction 
force . 

Impact  springs  were  incorporated  into  the  analytical  model  to  take 

account  of  joint  collisions.   The  stiffness  of  the  impact  spring  k 

7 
was  taken  equal  to  10  lbs/inch  which  is  nearly  the  same  magnitude  as 

the  longitudinal  stiffnesses  of  the  curved  beam  elements  representing 

the  superstructure.   The  initial  joint  gaps  A   were  determined  from 

G 

the  measured  experimental  responses  of  expansion  joints.   Since  the 
joint  gaps  were  changing  with  each  collision  and  also  were  not  identical 
between  the  inner  and  outer  edges,  the  following  joint  gaps  were 
obtained  by  taking  averages 


A  =  0.025  inch     expansion  joint  No.  1 
G  (84) 

A  =  0.040  inch     expansion  joint  No.  2 

Finally,  initial  tie  gaps  A   were  prescribed  on  the  basis  of 
measured  information  as 


A  =  0.01  inch     expansion  joint  No.  1 
T  (85) 

A  =  0.02  inch     expansion  joint  No.  2 

Since  no  experimental  data  was  available  to  estimate  damping 
factors  for  the  experimental  model  under  high  intensity  excitations, 
the  same  damping  coefficients  a  and  3  used  for  the  correlation  of 
Test  HI  were  adopted  as  a  first  trial.   As  before,  these  values 
represent  6%  viscous  damping  in  the  first  and  second  normal  modes  of 
vibration. 

After  introducing  these  data  into  the  nonlinear  analytical  model, 
it  was  subjected  to  the  table  accelerations  in  the  same  manner  as  done 
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previously  for  the  linear  analysis  correlation.   The  equilibrium  itera- 
tion was  used  both  alone  and  with  the  subdivision  of  time  intervals  into 
5  equidistant  sub-time  intervals,  when  needed.   The  tolerances  used  to 
control  the  equilibrium  iteration  and  the  subdivision  of  time  interval 
were 

A    =  0.01     equilibrium  iteration 
ps 

A    =  0.07     subdivision  (86) 

ps  v   ' 

A1   =  0.005 
us 

From  this  first  correlation,  it  was  found  that  although  the  overall 
response  was  substantially  improved  as  compared  with  the  linear  analysis 
correlation,  it  was  not  satisfactory.   The  discrepancies  were  primarily 
derived  from  excessive  constraint  of  the  tie  bars  at  both  expansion 
joints  and  from  over  estimation  of  the  damping  force.   The  former  reason 
indicated  that  the  assumed  initial  tie  gaps  were  too  small;  therefore, 
similar  to  the  correlation  of  Test  HI,  the  assumed  initial  tie  gaps 
were  increased.   In  order  to  interpret  the  latter  reason,  all  three 
sources  of  damping  in  the  analytical  model  were  examined,  i.e.  the  6% 
viscous  damping  which  was  proportional  to  both  mass  distribution  and 
stiffness  distribution,  Coulomb  type  damping  associated  with  slippages 
at  the  expansion  joints ,  and  hysteretic  damping  associated  with  yielding 
of  the  tie  bars.   The  6%  viscous  damping  was  originally  derived  from  the 
results  of  small  amplitude  forced  harmonic  excitations.   The  actual 
damping  in  these  tests  however  was  probably  more  in  the  form  of 
structural  damping  due  to  the  form  of  energy  dissipation  at  the  expansion 
joints  and  abutments.   Particularly  in  view  of  the  fact  that  damping  in 
the  sound  microconcrete  of  the  superstructure  and  columns  would  be  small. 
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Therefore,  when  collisions  take  place  and  high  frequency  vibration 
response  is  induced,  modelling  of  damping  using  the  same  6%  of  viscous 
damping  most  likely  would  overestimate  this  form  of  energy  dissipation. 

Thus,  changing  both  the  assumed  initial  tie  gaps  and  the  viscous 
damping  ratio,  a  best  correlation  was  achieved  under  the  conditions 

£  =  2%  of  critical  for  first  and  second  modes  of  vibration    (87) 

and 

A  =  0.03  inch     expansion  joint  No.  1 

(88) 

A  =  0.045  inch    expansion  joint  No.  2 

Comparative  plots  of  the  predicted  response  with  the  measured  response 
are  shown  in  Fig.  6.5.   Apparently  the  nonlinear  analysis  correlation  is 
greatly  superior  to  the  linear  analysis  correlation.   Clearly,  the  effect 
of  multiple  collisions  and  constraints  of  the  joint  restrainer  tie  bars 
are  realistically  represented  in  the  nonlinear  analytical  model,  see 
Fig.  6.6.   Severe  yielding  of  the  tie  bars  developed  maximum  ductility 
factors  u  ,   defined  as  the  ratio  of  the  maximum  tie  elongation  to  its 
yield  elongation,  i.e., 

E      P 

Um     +     Um 

^T  "      E  (89) 

UT 


where 

u   :   elastic  elongation  of  tie  bar 
T 

p 
u   :   plastic  elongation  of  tie  bar 

equal  to  approximately  4  and  12  for  expansion  joints  No.  1  and  No.  2, 
respectively.   Figure  6.7  shows  the  contact  force  between  girders  caused 
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by  collisions.   The  maximum  contact  force  of  approximately  4000  lbs  was 
induced  at  the  expansion  joint  No.  2.  Predicted  stresses  of  the  columns 
did  not  reach  to  their  yield  strength;  therefore,  no  plastic  deformation 
was  developed  during  the  test. 

C.   HIGH  INTENSITY  SEISMIC  EXCITATION  (TEST  H3) 

1.  Linear  Analytical  Correlation  -  Since  it  was  apparent  from  the 
correlation  of  Test  HV2  that  the  linear  analysis  could  not  achieve  good 
results  for  high  intensity  excitations,  the  linear  analysis  correlation 
was  performed  only  to  provide  comparative  response  data  with  the  non- 
linear analysis  correlation. 

The  same  linear  analytical  model  used  for  Test  HV2  was  used  for 
Test  H3.   It  was  subjected  to  the  table  motion  as  measured  in  the 
transverse  direction.   A  time  interval  of  0.01025  second  was  adopted 
throughout  this  correlation.   The  comparative  plots  of  the  predicted  and 
measured  response  are  shown  in  Fig.  6.8.   Like  test  HV2,  the  linearly 
predicted  response  is  very  much  different  from  the  measured  response. 

2.  Nonlinear  Analysis  Correlation  -  Since  the  general  behavior 

of  responses  in  this  excitation  test  is  essentially  the  same  as  in  Test 
HV2,  the  nonlinear  analytical  model  formulated  for  Test  HV2  was  again 
used.   The  same  damping  factors,  Coulomb  friction  coefficient,  impact 
springs  and  tie  bars  were  employed.   The  initial  tie  and  joint  gaps  were 
chosen  in  the  same  manner  described  previously  giving 

A  =  0.05  inch,    A  =  0.04  inch   expansion  joint  No.  1 
T  G  (90) 

A  =  0.1  inch,     A  =  0.05  inch   expansion  joint  No.  2 
T  G 

Comparative  plots  of  the  predicted  and  the  measured  responses  are  shown 
in  Fig.  6.9.   Like  Test  HV2 ,  the  nonlinear  analysis  correlation  is 
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greatly  superior  to  the  linear  analysis  correlation.   Although  permanent 
drifts  take  place  in  the  measured  response  from  approximately  2.1  seconds, 
it  is  beyond  the  scope  of  this  analytical  correlation  to  predict  them 
as  explained  in  Chapter  IV.   By  comparing  the  predicted  and  measured 
responses  at  expansion  joint  No.  2,  it  is  assessed  that  the  critical 
failure  of  the  shear  key  most  likely  was  caused  by  a  collision  which 
took  place  at  approximately  2  seconds.   Figure  6.10  shows  the  predicted 
tie  bar  forces  and  tie  yieldings  of  both  expansion  joints.   The  tie  bar 
ductility  factor  U   defined  in  Eq.   (89)   reached  approximately  10  and 
13  at  expansion  joints  No.  1  and  No.  2,  respectively.   Figure  6.11  shows 
the  contact  force  between  girders  caused  by  collisions.   A  maximum 
contact  force  of  approximately  8000  lbs  was  induced  at  expansion  joint 
No.  2.   Like  Test  HV2,  yielding  of  columns  did  not  develop  during  the 
test. 
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FIG.    6.1.1      LINEAR  CORRELATION   FOR  TEST   HI;    SIDE   GIRDER  NO.    1; 
LONGITUDINAL    (X)     DIRECTION 
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FIG.    6.1.2      LINEAR  CORRELATION   FOR  TEST  Hi;    SIDE   GIRDER  NO.    1; 
TRANSVERSE    (Y)    DIRECTION 
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FIG.  6.1.3   LINEAR  CORRELATION  FOR  TEST  HI;  CENTER  GIRDER; 
TRANSVERSE  (Y)  DIRECTION 
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FIG.  6.1.4   LINEAR  CORRELATION  FOR  TEST  HI;  SIDE  GIRDER  NO.  2; 
LONGITUDINAL  (X)  DIRECTION 
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FIG.  6.1.5   LINEAR  CORRELATION  FOR  TEST  HI;  SIDE  GIRDER  NO.  2; 
TRANSVERSE  (Y)  DIRECTION 


^0uZ3 


-  .  075 


1.0  1.5  2.0  2.5 

TIME      C  SECOND  ) 


3.0  3.5  4  .  0 


FIG.  6.1.6  LINEAR  CORRELATION  FOR  TEST  HI;  EXPANSION  JOINT  NO.  1; 
INNER  SIDE 
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FIG.  6.1.7   LINEAR  CORRELATION  FOR  TEST  Hi;  EXPANSION  JOINT  NO.  1; 
OUTER  SIDE 
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FIG.  6.1.8   LINEAR  CORRELATION  FOR  TEST  HI;  EXPANSION  JOINT  NO.  2; 
INNER  SIDE 
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FIG.  6.1.9   LINEAR  CORRELATION  FOR  TEST  Hi;  EXPANSION  JOINT  NO.  2; 
OUTER  SIDE 
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FIG.  6.2.1   NONLINEAR  CORRELATION  FOR  TEST  HI;  SIDE  GIRDER  NO.  1; 
LONGITUDINAL  (X)  DIRECTION 
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FIG.  6.2.2   NONLINEAR  CORRELATION  FOR  TEST  Hi;  SIDE  GIRDER  NO.  1; 
TRANSVERSE  (Y)  DIRECTION 
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FIG.  6.2.3   NONLINEAR  CORRELATION  FOR  TEST  HI;  CENTER  GIRDER; 
TRANSVERSE  (Y)  DIRECTION 


.07  5 


1.0        1.5        2.0        2.5 

TIME   (SECOND] 


3.0        3.5        «  .  0 


FIG.  6.2.4   NONLINEAR  CORRELATION  FOR  TEST  HI;  SIDE  GIRDER  NO.  2; 
LONGITUDINAL  (X)  DIRECTION 
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FIG.  6.2.5   NONLINEAR  CORRELATION  FOR  TEST  Hi;  SIDE  GIRDER  NO.  2; 
TRANSVERSE  (Y)  DIRECTION 
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FIG.  6.2.6   NONLINEAR  CORRELATION  FOR  TEST  HI;  EXPANSION  JOINT  NO.  1; 
INNER  SIDE 
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FIG.  6.2.7   NONLINEAR  CORRELATION  FOR  TEST  HI;  EXPANSION  JOINT  NO.  1; 
OUTER  SIDE 
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FIG.  6.2.8   NONLINEAR  CORRELATION  FOR  TEST  HI;  EXPANSION  JOINT  NO.  2; 
INNER  SIDE 
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FIG.  6.2.9   NONLINEAR  CORRELATION  FOR  TEST  HI;  EXPANSION  JOINT  NO.  2; 
OUTER  SIDE 
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FIG.  6.3.1   PREDICTED  TIE  BAR  FORCE  FOR  TEST  Hi;  EXPANSION  JOINT  NO.  1 
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FIG.  6.3.2   PREDICTED  TIE  BAR  FORCE  FOR  TEST  HI;  EXPANSION  JOINT  NO.  2 
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FIG.  6.4.1   LINEAR  CORRELATION  FOR  TEST  HV2 ;  SIDE  GIRDER  NO.  1; 
LONGITUDINAL  (X)  DIRECTION 
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FIG.  6.4.2   LINEAR  CORRELATION  FOR  TEST  HV2 ;  SIDE  GIRDER  NO.  1  ; 
TRANSVERSE  (Y)  DIRECTION 
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FIG.  6.4.3   LINEAR  CORRELATION  FOR  TEST  HV2 ;  CENTER  GIRDER; 
TRANSVERSE  (Y)  DIRECTION 
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FIG.  6.4.4   LINEAR  CORRELATION  FOR  TEST  HV2 ;  SIDE  GIRDER  NO.  2 ; 
LONGITUDINAL  (X)  DIRECTION 
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FIG.  6.4.5   LINEAR  CORRELATION  FOR  TEST  HV2 ;  SIDE  GIRDER  NO.  2 ; 
TRANSVERSE  (Y)  DIRECTION 
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FIG.  6.4.6   LINEAR  CORRELATION  FOR  TEST  HV2 ;  EXPANSION  JOINT  NO.  1; 
INNER  SIDE 
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FIG.  6.4.7   LINEAR  CORRELATION  FOR  TEST  HV2 ;  EXPANSION  JOINT  NO.  1 ; 
OUTER  SIDE 
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FIG.  6.4.8   LINEAR  CORRELATION  FOR  TEST  HV2  -r    EXPANSION  JOINT  NO.  2; 
INNER  SIDE 
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FIG.  6.4.9   LINEAR  CORRELATION  FOR  TEST  HV2 ;  EXPANSION  JOINT  NO.  2 ; 
OUTER  SIDE 
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FIG.  6.5.1   NONLINEAR  CORRELATION  FOR  TEST  HV2 ;  SIDE  GIRDER  NO.  1 
LONGITUDINAL  (X)  DIRECTION 
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FIG.  6.5.2   NONLINEAR  CORRELATION  FOR  TEST  HV2 ;  SIDE  GIRDER  NO.  1 ; 
TRANSVERSE  (Y)  DIRECTION 

117 


3  50 


2  33 


h- 


UJ    .117 

E 


Q. 
C/) 


.117 


<  -.233 


3  5  0 


Predicted 
Measured 


.5  1.0 


1.5         2.0         2.5 

TIME   (SECOND) 


3.0         3.5         4.0 


FIG.  6.5.3   NONLINEAR  CORRELATION  FOR  TEST  HV2 ;  CENTER  GIRDER; 
TRANSVERSE  (Y)  DIRECTION 
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FIG.  6.5.4   NONLINEAR  COORELATION  FOR  TEST  HV2 ;  SIDE  GIRDER  NO.  2; 
LONGITUDINAL  (X)  DIRECTION 
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FIG.  6.5.5   NONLINEAR  CORRELATION  FOR  TEST  HV2 ;  SIDE  GIRDER  NO.  2  ; 
TRANSVERSE  (Y)  DIRECTION 
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FIG.  6.5.6   NONLINEAR  CORRELATION  FOR  TEST  HV2 ;  EXPANSION  JOINT  NO.  1; 
INNER  SIDE 
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FIG.  6.5.7   NONLINEAR  CORRELATION  FOR  TEST  HV2 ;  EXPANSION  JOINT  NO.  1; 
OUTER  SIDE 
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FIG.  6.5.8   NONLINEAR  CORRELATION  FOR  TEST  HV2 ;  EXPANSION  JOINT  NO.  2; 
INNER  SIDE 
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FIG.  6.5.9   NONLINEAR  CORRELATION  FOR  TEST  HV2 ;  EXPANSION  JOINT  NO.  2; 
OUTER  SIDE 
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FIG.  6.6.1   PREDICTED  TIE  BAR  FORCE  AND  YIELD  LENGTH  FOR  TEST  HV2 ; 
EXPANSION  JOINT  NO.  1 
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FIG.  6.6.2   PREDICTED  TIE  BAR  FORCE  AND  YIELD  LENGTH  FOR  TEST  HV2 ; 
EXPANSION  JOINT  NO.  2 
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FIG.  6.7.1   PREDICTED  CONTACT  FORCE  FOR  TEST  HV2 ;  EXPANSION  JOINT  NO.  1 
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FIG.  6.7.2   PREDICTED  CONTACT  FORCE  FOR  TEST  HV2 ;  EXPANSION  JOINT  NO.  2 
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FIG.  6.8.1   LINEAR  CORRELATION  FOR  TEST  H3;  SIDE  GIRDER  NO.  1; 
LONGITUDINAL  (X)  DIRECTION 
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FIG.  6.8.2   LINEAR  CORRELATION  FOR  TEST  H3;  SIDE  GIRDER  NO.  1; 
TRANSVERSE  (Y)  DIRECTION 
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FIG.  6.8.3   LINEAR  CORRELATION  FOR  TEST  H3;  CENTER  GIRDER; 
TRANSVERSE  (Y)  DIRECTION 


.4  00 


c_>        .300 


r 

UJ 

< 

Q- 


.20  0 


100     ■ 


100     ■ 


2  00     ■ 


-  .  300 


4  00 


1  .  0 


1  .  5 

TIME 


2.0  2.5 

[SECOND) 


3  .  0 


3  .  5 


4  .  0 


FIG.  6.8.4   LINEAR  CORRELATION  FOR  TEST  H3;  SIDE  GIRDER  NO.  2 ; 
LONGITUDINAL  (X)  DIRECTION 
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FIG.  6.8.5   LINEAR  CORRELATION  FOR  TEST  H3;  SIDE  GIRDER  NO.  2 ; 
TRANSVERSE  (Y)  DIRECTION 
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FIG.  6.8.6   LINEAR  CORRELATION  FOR  TEST  H3;  EXPANSION  JOINT  NO.  1; 
INNER  SIDE 
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FIG.  6.8.7   LINEAR  CORRELATION  FOR  TEST  H3;  EXPANSION  JOINT  NO.  1 ; 
OUTER  SIDE 
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FIG.  6.8.8   LINEAR  CORRELATION  FOR  TEST  H3 ;  EXPANSION  JOINT  NO.  2; 
INNER  SIDE 
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FIG.  6.8.9   LINEAR  CORRELATION  FOR  TEST  H3 ;  EXPANSION  JOINT  NO.  2; 
OUTER  SIDE 
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FIG.  6.9.1   NONLINEAR  CORRELATION  FOR  TEST  H3;  SIDE  GIRDER  NO.  1; 
LONGITUDINAL  (X)  DIRECTION 


100 


-  .  100 


1  .  0 


1  .  5 

TIME 


2.0        2.5 

(SECOND] 


3  .  o 


3  .  5 


1  .  0 


FIG.  6.9.2   NONLINEAR  CORRELATION  FOR  TEST  H3 ;  SIDE  GIRDER  NO.  1; 
TRANSVERSE  (Y)  DIRECTION 
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FIG.  6.9.3   NONLINEAR  CORRELATION  FOR  TEST  H3 ;  CENTER  GIRDER; 
TRANSVERSE  (Y)  DIRECTION 
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FIG.  6.9.4   NONLINEAR  CORRELATION  FOR  TEST  H3 ;  SIDE  GIRDER  NO.  2 ; 
LONGITUDINAL  (X)  DIRECTION 
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FIG.  6.9.5   NONLINEAR  CORRELATION  FOR  TEST  H3 ;  SIDE  GIRDER  NO.  2 ; 
TRANSVERSE  (Y)  DIRECTION 
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FIG.  6.9.6   NONLINEAR  CORRELATION  FOR  TEST  H3;  EXPANSION  JOINT  NO.  1; 
INNER  SIDE 
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FIG.  6.9.7   NONLINEAR  CORRELATION  FOR  TEST  H3 ;  EXPANSION  JOINT  NO.  1; 
OUTER  SIDE 


lj   .3  00 


.400 


1.5        2.0        2.5 

TIME       (SECOND] 


Predicted 
Measured 


3  .  o 


FIG.  6.9.8   NONLINEAR  CORRELATION  FOR  TEST  H3 ;  EXPANSION  JOINT  NO.  2 ; 
INNER  SIDE 


132 


4  0  0 


1  .  0 


1.5         2.0         2.5 

TIME   (SECOND] 


3  .  0 


3  .  5 


4  .  0 


FIG.  6.9.9   NONLINEAR  CORRELATION  FOR  TEST  H3;  EXPANSION  JOINT  NO.  2; 
OUTER  SIDE 
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FIG.  6.10.1   PREDICTED  TIE  BAR  FORCE  AND  YIELD  LENGTH  FOR  TEST  H3; 
EXPANSION  JOINT  NO.  1 
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VII.   DISCUSSION  OF  CORRELATION  RESULTS 

The  purpose  of  this  chapter  is  to  comment  on  the  correlation  results 
presented  in  the  previous  chapter.   Since  Test  HV2  shows  the  most 
typical  response  of  the  bridge  model  under  high  intensity  excitation 
including  the  tie  bar  yieldings, collisions,  slippages  and  vertical 
excitation,  these  comments  will  concentrate  on  this  test. 

A.   JOINT  RESTRAINER  TIE  BARS 

The  experimental  model  consisted  of  three  independent  subassemblages 
having  essentially  different  dynamic  characteristics  which  were  joined 
together  by  restrainer  tie  bars.   The  characteristics  of  these  bars  are 
extremely  important  in  controlling  the  dynamic  behavior  of  the  bridge 
model  both  in  the  amplitude  and  in  frequency  content.   Figure  7.1  shows 
an  example  of  two  different  calculated  responses  which  are  compared  with 
the  corresponding  predicted  response  adopted  for  the  correlation  of 
Test  HV2.   These  three  responses  were  calculated  using  the  same  nonlinear 
analytical  model  and  damping  factors  but  with  different  tie  bar  char- 
acteristics.  In  one  case,  the  magnitude  of  the  initial  tie  gap  at 
expansion  joint  No.  2  was  arbitrarily  taken  twice  that  value  adopted  for 
the  correlation.   In  another  case,  it  was  assumed  that  no  tie  bars  were 
provided  at  the  expansion  joints.   It  is  apparent  from  this  comparison 
that  the  joint  restrainer  tie  bars  can  be  extremely  effective  in 
reducing  the  seismic  response  of  the  bridge  model,  if  properly  designed. 
Further,  it  is  apparent  that  small  initial  tie  gaps  are  desirable  for 
this  purpose. 
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B.  COLLISIONS 

Although  the  general  effects  of  collisions  on  the  dynamic  behavior 
of  the  bridge  model  is  apparent  from  a  comparison  of  responses  between 
the  linear  analysis  correlation  and  the  nonlinear  analysis  correlation 
for  the  Test  HV2  or  Test  H3,  the  differences  observed  in  this  comparison 
include  the  effects  of  other  factors,  e.g.  tie  bars,  slippages,  etc. 
Therefore  in  an  attempt  to  isolate  the  effects  of  collisions,  one  response 
was  calculated  using  the  same  nonlinear  analytical  model  adopted  finally 
for  the  correlation  of  the  Test  HV2  but  eliminating  the  impact  springs 
at  both  expansion  joints.   The  comparative  plot  between  this  response 
and  the  response  calculated  for  the  correlation  of  Test  HV2  is  shown  in 
Fig.  7.2.   It  is  seen  from  this  comparison  that  collisions  are  a  major 
factor  affecting  the  general  behavior  of  the  bridge  model  subjected  to 
high  intensity  seismic  excitations.   It  should  be  noted  that  even  if  the 
response  amplitude  in  the  negative  direction  is  suppressed  by  a  collision, 
the  rebound  from  that  collision  can  appreciably  increase  the  positive 
response  amplitude  which  follows  immediately. 

C.  COULOMB  FRICTION  FORCE 

Although  a  Coulomb  friction  coefficient  of  0.4  was  assumed  through- 
out this  study  based  on  considerations  of  the  surface  condition  of  the 
contact  plane,  its  actual  value  could  be  expected  to  vary  from  0.3  to  0.6. 
Therefore,  it  is  of  interest  to  examine  the  sensitivity  of  this  para- 
meter on  the  dynamic  behavior  of  the  bridge  model.   Therefore,  two 
responses  were  calculated  using  the  nonlinear  analytical  model  formulated 
for  Test  HV2  but  with  different  Coulomb  friction  coefficients,  i.e.  0.3 
and  0.6.   The  comparative  plot  of  these  two  responses  is  shown  in 
Fig.  7.3  which  could  be  compared  with  the  response  finally  adopted  in 
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the  correlation  of  Test  HV2  using  a  friction  coefficient  of  0.4.   It  is 
apparent  from  this  comparison  that  the  Coulomb  friction  coefficient  has 
low  sensitivity  in  controlling  seismic  response.   For  example,  if  a 
Coulomb  friction  coefficient  of  0.6  is  used  which  could  be  considered 
a  maximum,  the  maximum  friction  force  acting  on  one  expansion  joint 
would  be  only  220  lbs  corresponding  to  a  vertical  contact  force  of 
380  lbs.   In  comparison,  the  maximum  tie  bar  force  acting  on  one  expan- 
sion joint  is  approximately  1300  lbs  which  is  roughly  6  times  as  large. 
Furthermore,  the  predicted  maximum  contact  force  caused  by  collisions 
reached  approximately  4000  lbs  which  is  roughly  18  times  as  large.   It 
is  clear  from  this  comparison  that  the  overall  response  was  critically 
controlled  by  the  tie  bars  for  outward  motion  and  by  collisions  for 
inward  motion  with  the  effects  of  friction  being  relatively  small. 

D.   VERTICAL  EXCITATION 

Vertical  excitation  has  an  important  effect  on  the  vertical  res- 
ponse of  the  bridge  model.   It  also  has  some  effect  on  horizontal  res- 
ponse due  to  coupling  between  the  horizontal  and  vertical  vibration 
modes  which  is  associated  with  deck  curvature.   Further,  the  vertical 
excitation  causes  vertical  oscilatory  motion  in  the  superstructure 
resulting  in  changes  of  the  vertical  contact  force  at  the  expansion 
joints.   This,  in  turn,  affects  the  friction  forces  acting  at  expansion 
joints  which  affect  horizontal  response. 

To  check  the  influence  of  vertical  excitation,  one  response  was 
calculated  using  the  nonlinear  analytical  model  adopted  for  the  cor- 
relation of  Test  HV2  but  with  the  vertical  excitation  removed  from  the 
input.   The  comparative  plot  of  responses  with  and  without  vertical 
excitation  is  shown  in  Fig.  7.4.   It  is  apparent  from  this  comparison 
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that  vertical  excitation  is  less  significant  than  horizontal  excitation 
in  the  transverse  direction.   This  observation  is  consistent  with  the 
relatively  low  sensitivity  of  Coulomb  friction  and  relatively  small 
effect  caused  by  coupling  of  vertical  and  horizontal  modes. 

E.   IMPACT  SPRING  STIFFNESS 

In  the  analytical  correlations  presented  in  Chapter  VI,  the  impact 

7 
spring  stiffness   k   of  10   lbs/inch  was  adopted.   This  value  was 

selected  so  that  k   would  be  similar  in  magnitude  to  the  longitudinal 
stiffness  of  neighboring  curved  beam  elements  representing  the  super- 
structure; see  Appendix  A.   Since  impact  is  important  to  response,  the 
sensitivity  of  stiffness  k   was  examined  using  the  nonlinear  analytical 

model  adopted  for  the  correlation  of  Test  HV2  but  using  two  different 

6       9 

impact  spring  stiffnesses  k  ,   namely  10  and  10   lbs/inch  which  are, 

respectively,  one  tenth  and  a  hundred  times  that  of  the  impact  spring 
stiffness  adopted  for  the  correlation  study.   It  is  apparent  from  the 
resulting  responses  shown  in  Fig.  7.5  that  k   is  an  insensitive  para- 
meter provided  it  is  set  sufficiently  large.   Therefore,  the  stiffness 

7 
of  10   lbs/inch  set  in  the  correlation  study  is  shown  to  be  a  proper 

value  for  the  analysis  purpose. 
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VIII   CONCLUSIONS 

Based  on  the  correlation  results  presented,  the  following  conclusions 
may  be  deduced: 

(1)  The  mathematical  model  and  analytical  procedure  developed  in 
Phase  2 ,  which  were  modified  and  improved  in  Phase  5 ,  predicts  realis- 
tically the  nonlinear  response  of  the  bridge  model  under  low  and  high 
intensity  seismic  excitations  as  conducted  in  Phase  4. 

(2)  The  displacement  seismic  response  of  the  bridge  model  under  low 
intensity  excitation  which  does  not  result  in  significant  yielding  of 
the  longitudinal  joint  restrainer  tie  bars  and  does  not  produce  collision 
of  the  girders  can  be  predicted  with  fairly  good  accuracy  using  the 
linear  analytical  model  provided  the  low  frequency  characteristics  of 
analytical  model  agree  well  with  those  of  the  physical  model  and  pro- 
vided a  proper  damping  factor  is  used  in  the  analysis. 

(3)  The  linear  analytical  model  cannot  satis f actor ly  predict  the 
seismic  response  of  the  bridge  model  under  high  intensity  seismic 
excitations  causing  appreciable  yielding  of  tie  bars  and  multiple 
collisions  of  the  girders,  even  if  the  low  frequency  characteristics 
and  the  damping  factor  of  the  analytical  model  are  adjusted. 

(4)  The  displacement  seismic  response  of  the  bridge  model  produced 
by  high  intensity  excitations  resulting  in  severe  yielding  of  tie  bar 
and  collisions  of  the  girders  can  be  predicted  realistically  using  the 
nonlinear  analytical  model  which  accounts  for  the  effects  of  collisions, 
slippages  and  tie  bars. 

(5)  The  expansion  joints  have  a  controlling  effect  on  the  dynamic 
behavior  of  the  bridge  model;  thus,  their  characteristics  must  be  pro- 
perly represented  in  a  realistic  ma thematic  model  of  the  bridge  structure. 
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(6)  Under  high  intensity  seismic  excitations,  the  dynamic  response 
of  the  bridge  model  is  primarily  controlled  by  joint  restrainer  tie  bars 
for  motions  in  the  outward  direction  and  by  multiple  collisions  for 
motions  in  the  inward  direction  with  the  effects  of  Coulomb  friction 
being  relatively  small. 

(7)  Vertical  excitation  has  an  important  effect  on  the  vertical 
response  of  the  bridge  model.   However,  the  effect  of  vertical  excitation 
on  horizontal  transverse  response  is  relatively  small  due  to  the 
insensitivity  of  coupling  of  vertical  and  horizontal  modes  and  of  Coulomb 
friction  on  horizontal  response. 

(8)  The  longitudinal  restrainer  tie  bars  are  very  effective  in 
reducing  transverse  response  of  the  bridge  model  provided  they  are  pro- 
perly designed.   Smaller  tie  gaps  are  also  effective  in  reducing  trans- 
verse response. 

(9)  The  multiple  collisions  which  take  place  between  girders  have 

a  major  influence  on  both  the  amplitude  and  the  frequency  characteristics 
of  response  and  they  cause  large  contact  forces  to  be  developed  at 
expansion  joints.   Even  though  a  negative  response  amplitude  is  sup- 
pressed by  collision,  the  rebound  from  that  collision  can  result  in  an 
increased  positive  response  amplitude  immediately  following  the  collision. 

(10)  The  mathematical  model,  analytical  procedures,  and  computer 
program  NEABS  can  be  used  to  predict  realistic  dynamic  response  of  proto- 
type bridge  structures  subjected  to  high  intensity  rigid  base  seismic 
excitation. 


146 


BIBLIOGRAPHY 


1.  Iwasaki,  T. ,  Penzien,  J.  and  Clough,  R. ,  "Literature  Survey  — 
Seismic  Effects  on  Highway  Bridges,"  Report  No.  EERC  71-11,  Earth- 
quake Engineering  Research  Center,  University  of  California, 
Berkeley,  November  1972. 

2.  Tseng,  W.  S.,  and  Penzien,  J.,  "Analytical  Investigations  of  the 
Seismic  Response  of  Long  Multiple  Span  Highway  Bridges , "  Report 
No.  EERC  73-12,  Earthquake  Engineering  Research  Center,  University 
of  California,  Berkeley,  June  1973. 

3.  Tseng,  W.  S.,  and  Penzien,  J.,  "Seismic  Response  of  Highway  Over- 
crossings,"  Proc.  5th  World  Conference  on  Earthquake  Engineering, 
Rome,  Italy,  1973. 

4.  Tseng,  W.  S.,  and  Penzien,  J.,  "Seismic  Analysis  of  Long  Multiple- 
Span  Highway  Bridges,"  International  Journal  of  Earthquake  Engineering 
and  Structural  Dynamics,  Vol.  4,  pp.  3-24,  1975. 

5.  Tseng,  W.  S.,  and  Penzien,  J.,  "Seismic  Response  of  Long  Multiple- 
Span  Highway  Bridges,"  International  Journal  of  Earthquake 
Engineering  and  Structural  Dynamics,  Vol.  4,  pp.  25-48,  1975. 

6.  Chen,  M.  C. ,  and  Penzien,  J.,  "Analytical  Investigations  of  Seismic 
Response  of  Short,  Single,  or  Multiple-Span  Highway  Bridges,^  Report 
No.  EERC  7  5-4,  Earthquake  Engineering  Research  Center,  University  of 
California,  Berkeley,  January  1975. 

7.  Goldsmith,  W. ,  "Impacts,"  Edward  Arnold,  London,  1960. 

8.  Clough,  R.  W. ,  and  Penzien,  J.,  "Dynamics  of  Structure,  McGraw  Hill, 
1975. 

9.  Ban,  S.,  "Collisions  of  Building  Structure  during  Earthquakes," 
Trans,  of  the  Architectural  Institute  of  Japan,  Vol.  221,  pp.  1-7, 
July,  1974,  (In  Japanese). 

10.  Takeyama,  K. ,  "Earthquake  Response  of  A  Building  Collied  with  A 
Neighboring  Building,"  Proc.  of  5th  World  Conference  on  Earthquake 
Engineering,  Rome,  Italy,  197  3. 

11.  Bathe,  K.  J.,  Wilson,  E.  L. ,  and  Peterson,  F.  E.,  "SAP  IV  —  A 
Structural  Analysis  Program  for  Static  and  Dynamic  Response  of 
Linear  Systems,  Report  No.  EERC  73-11,  Earthquake  Engineering 
Research  Center,  University  of  California,  Berkeley,  June  1973. 

12.  Newmark,  N.  M. ,  "A  Method  of  Computation  for  Structural  Dynamics," 
Proc.  ASCE,  Vol.  85,  No.  EM3 ,  pp.  67-94,  July  1959. 

13.  Wilson,  E.  L.  ,  and  Clough,  R.  W.  ,  ''Dynamic  Response  by  Step-by-Step 
Matrix  Analysis,"  Proc.  Symposium  on  the  use  of  Computers  in  Civil 
Engineering,  Lisbon,  Portugal,  1962. 

147 


14.  Clough,  R.  W. ,  and  Bathe,  K.  J.,  "Finite  Element  Analysis  of  Dynamic 
Response,"  2nd  U.S. -Japan  Seminar  on  Matrix  Methods  of  Structural 
Analysis  and  Design,  pp.  153-179,  August  1972. 

15.  Bathe,  K.  J.,  and  Wilson,  E.  L. ,  "Stability  and  Accuracy  Analysis 
of  Direct  Integration  Methods,"  International  Journal  of  Earthquake 
Engineering  and  Structural  Dynamics,  Vol.  1,  pp.  282-291,  1973. 

16.  Wilson,  E.  L. ,  Farhoomand,  I.,  and  Bathe,  K.  J.,  "Nonlinear  Dynamic 
Analysis  of  Complex  Structures,"  International  Journal  of  Earthquake 
Engineering  and  Structural  Dynamics,  Vol.  1,  pp.  241-252,  1973. 

17.  Nickell,  R.  E.,  "Direct  Integration  Methods  in  Structural  Dynamics," 
Proc.  ASCE,  Vol.  99,  No.  EM2,  pp.  303-317,  April  1973. 

18.  Kanaan,  A.  E. ,  and  Powell,  G.  H. ,  "General  Purpose  Computer  Program 
for  Dynamics  Analysis  of  Inelastic  Plane  Structures ,"  Report  No.  EERC 
73-6,  Earthquake  Engineering  Research  Center,  University  of 
California,  Berkeley,  1973. 

19.  Hughes,  T.  J.,  Taylor,  R.  L. ,  and  Sackman,  J.  L. ,  "Finite  Element 
Formulation  and  Solution  of  Contact  -  Impact  Problems  in  Continuous 
Mechanics,"  Report  No.  UC  SESM  74-8,  Structural  Engineering 
Laboratory,  University  of  California,  Berkeley,  1974. 

20.  Mondkar,  D.  P.,  and  Powell,  G.  H. ,  "Static  and  Dynamic  Analysis  of 
Nonlinear  Structures,"  Report  No.  EERC  75-10,  University  of 
California,  Berkeley,  March  1975. 

21.  Bathe,  K.  J.,  Wilson,  E.  L. ,  and  Iding,  R.  H. ,  "NONSAP  —  A 
Structural  Analysis  Program  for  Static  and  Dynamic  Response  of 
Nonlinear  System,"  Report  No.  UC  SESM  74-3,  Structural  Engineering 
Laboratory,  University  of  California,  Berkeley,  1974. 

22.  Bathe,  K.  J.,  Ozdemir,  H. ,  and  Wilson,  E.  L. ,  "Static  and  Dynamic 
Geometric  and  Material  Nonlinear  Analysis,"  Report  No.  UC  SESM  74-4, 
Structural  Engineering  Laboratory,  University  of  California, 
Berkeley,  1974. 

23.  Williams,  D. ,  and  Godden,  W.  G. ,  "Seismic  Response  of  A  Curved 
Highway  Bridge  Model,"  54th  Annual  Meeting,  Transportation  Research 
Board,  National  Research  Council,  Washington,  D.  C. ,  January,  1975. 

24.  Williams,  D. ,  and  Godden,  W.  G. ,  "Seismic  Behavior  of  High  Curved 
Overcrossings , "  Proc.  of  U.  S.  National  Conference  on  Earthquake 
Engineering,  Ann  Arbor,  1975. 

25.  Howe,  S. ,  and  Williams,  D. ,  "Private  Memorandom  on  Test  Procedures 
and  Test  Results  of  Model  Bridge  Structure,"  1975. 

26.  Williams,  D. ,  and  Godden,  W.  G. ,  "16  m/m  Film  for  Experimental 
Test  of  Model  Bridge  Structure,  1975. 


148 


27.  Rea,  D. ,  and  Penzien,  J.,  "Dynamic  Response  of  A  20  ft  x  20  ft 
Shaking  Table,"  Proc.  of  5th  World  Conference  on  Earthquake 
Engineering,  Rome,  Italy,  1973. 

28.  Chopra,  A.  K. ,  Bertero,  V.  V.,  and  Mahin,  S. ,  "Response  of  the  Olive 
View  Medical  Center  Main  Building  during  the  San  Fernando  Earthquake," 
Proc.  5th  World  Conference  on  Earthquake  Engineering,  Rome,  Italy, 
1973. 

29.  Tang,  D. ,  and  Clough,  R.  W. ,  "Shaking  Table  Tests  of  A  Steel  Frame  - 
A  Progress  Report,"  Report  No.  EERC  74-8,  Earthquake  Engineering 
Research  Center,  University  of  California,  Berkeley,  July  1975. 

30.  Tang,  D. ,  "Earthquake  Simulatory  Study  of  A  Steel  Frame  Structure, 
Volume  II  -  Analytical  Results,"  Report  No.  EERC  75-36,  Earthquake 
Engineering  Research  Center,  University  of  California,  Berkeley, 
October  1975. 

31.  Hidalgo,  P.,  and  Clough,  R.  W. ,  "Earthquake  Simulatory  Study  of  A 
Reinforced  Concrete  Frame,"  Report  No.  EERC  74-13,  Earthquake 
Engineering  Research  Center,  University  of  California,  Berkeley, 
December  1974. 

32.  Iemura,  H. ,  and  Jennings,  P.  C. ,  "Hysteretic  Response  of  a  Nine-Story 
Reinforced  Concrete  Building  During  the  San  Fernando  Earthquake , " 
Report  No.  EERL  73-7,  Earthquake  Engineering  Research  Laboratory, 
California  Institute  of  Technology,  Pasadena,  October  1973. 

33.  Bouwkamp,  J.  G. ,  and  Rea,  D. ,  "Dynamic  Testing  and  the  Formulations 
of  Mathematical  Models,"  Earthquake  Engineering,  Wiegel ,  R.  L. , 
Coordinating  Editor,  Prentice  Hall,  1970. 


149 


APPENDIX  A 
LONGITUDINAL  COLLINEAR  COLLISION  OF  TWO  RODS 

To  check  the  mathematical  model  and  analytical  procedures 
used  in  representing  joint  impact  as  presented  in  Chapters  II  and  III, 
respectively,  consider  the  collinear  collision  of  two  uniform  bars 
initially  travelling  in  opposite  directions  with  the  same  initial 
velocity  V   as  shown  in  Fig.  A.l.   The  approximate  solution  of  the 
post- impact  behavior  obtained  by  these  procedures  can  be  compared  with 
the  exact  solution  obtained  by  classical  wave  propagation  theory  [7] . 

Assume  both  bars  have  the  same  properties  as  given  by 

E  (modulus  of  elasticity)  =  100 

A  (cross-sectional  area)  =  1 

p  (mass  density)  =  0.1  (A-l) 

L  (rod  length)  =10 

V   (initial  rod  velocity)  =  ±  0.1 
o 

in  which  any  convenient  units  may  be  used.  For  this  example  problem, 
the  impact  contact  duration  T  is  0.2  units  of  time  as  given  by  the 
exact  solution.  This  duration  corresponds  to  the  time  required  for  a 
wave  to  propagate  twice  the  length  of  the  rod,  i.e. 

T   =   —  (A-2) 

I     c 


where  c   is  the  longitudinal  wave  velocity  given  by 


c   =   /E/p  (A-3) 

o 

In  both  the  exact  and  approximate  solutions ,  it  is  convienient 
to  monitor  the  separation  and  rate  of  change  of  separation  between  the 
two  bars  as  given  by 


A-l 


R     L 
u  =  u  -  u 

-    -R    -L  (A'4) 

u  =  u  -  u 

R       *R 
where  u   and  u   represent  the  displacement  and  velocity,  respectively, 

L       »L 
of  the  contact  surface  of  the  right  bar  and  u   and  u   represent  the 

displacement  and  velocity,  respectively,  of  the  contact  surface  of  the 
left  bar.   A  positive  value  of  u  represents  a  separation  of  the  con- 
tact surfaces  while  a  negative  value  of  u  represents  an  overlap  of  the 
contact  surfaces.   A  negative  value  of  u  is,  of  course,  not  possible 
in  reality  as  given  by  the  exact  solution. 

Using  the  modelling  procedures  of  Chapter  II,  the  two  bars  can 
be  represented  by  an  impact  spring  and  finite  elements  as  shown  in 
Fig.  A-2.   In  this  case,  a  zero  value  of  u  represents  that  position  of 
the  bars  when  initial  contact  is  made  on  both  sides  with  the  impact 
spring;  thus,  a  negative  value  of  u  is  possible  which  corresponds  to 
the  shortening  of  the  impact  spring. 

In  modelling  the  two  bar  system  in  Fig.  A-2,  the  impact  spring 
stiffness  K  ,   the  numerical  time  interval  of  integration  At,   and  the 
number  of  finite  elements  n  representing  each  rod  are  all  factors 
influencing  the  predicted  dynamic  response  following  impact.   The  con- 
stant acceleration  integration  procedure  with  equilibrium  iteration  has 
been  used  in  this  example  parameter  study.   The  tolerances  controlling 
the  equilibrium  iteration,  as  defined  in  Eqs.  (51)  and  (63),  were 
specified  as 

A    =  0.001 
ps 

(A-5) 

A1   =  0.001 

ps 


A-2 


Responses  u  and  u  as  predicted  by  exact  wave  propagation 
theory  are  shown  in  Fig.  A. 3.   Before  collision,  the  relative  separation 
velocity  u  equals  -0.2  which  instantaneously  changes  to  zero  upon 
collision.   This  velocity  stays  at  a  value  of  zero  during  the  contact 
duration  equal  to  0.2  and  then  instantaneously  changes  to  a  value  of 
+0.2.   It  should  be  noted  that  the  instantaneous  change  in  velocity 
represents  a  Dirac  delta  change  in  the  acceleration,  i.e.  a  pure 
acceleration  pulse  of  duration  and  amplitude  which  tend  to  zero  and 
infinity,  respectively. 

Turning  now  to  the  approximate  solution,  consider  first  the 
effect  of  the  magnitude  of  the  impact  spring  stiffness  K  upon  res- 
ponse.  Using  a  time  interval  of  integration  At  equal  to  0.01  T  ,  a 

I        2    3 
value  of  10  for  n,   and  five  different  values  for  K   (10,  10  ,  10  , 

4       5 
10  and  10  ) ,  responses  u  and  u  as  shown  in  Fig.  A. 4  were  obtained. 

Comparing  the  relative  displacement  results  in  this  figure  with  the 

corresponding  results  for  the  exact  solution  given  in  Fig.  A. 3  shows 

that  the  exact  and  approximate  values  of  u  are  in  good  agreement  for 

I  2       3 

K   equal  to  10  and  10  but  they  are  significantly  in  error  for  the 

larger  values  of  K  .   Comparing  the  relative  velocity  results  in  these 

same  figures  shows  considerable  differences  for  all  5  values  of  K   and 

it  shows  considerable  vibrations  taking  place  following  unit  time  zero. 

These  vibrations  in  velocity  u  increase  in  amplitude  with  increasing 

values  of  K  .   It  is  interesting  to  note  that  when  the  larger  values 

of  K   are  used,  contact  and  separation  of  the  bar  ends  are  repeated 

many  times  during  period  T  .   To  judge  the  best  value  of  K   it  is 

convenient  to  introduce  the  dimensionless  parameter 

v  =  £__L  (A_6) 

1     n  EA 

A-3 


which  represents  the  ratio  of  the  impact  spring  stiffness  to  an 

individual  rod  element  stiffness .   Comparing  responses  u  and  u  in 

Figs.  A. 3  and  A. 4,  it  appears  that  a  value  of  unity  for  y     yields  the 

best  overall  correlation  of  the  approximate  results  with  exact  results. 

This  suggests  that  the  numerical  value  of  K   should  be  approximately 

equal  to  the  stiffness  of  its  neighboring  elements. 

Let  us  now  consider  the  effect  of  the  magnitude  of  At  upon 

response.   Using  K   equal  to  100  (7=1),   n  equal  to  10,  and  At 

equal  to  0.5  T  ,  0.05  T  ,  and  0.01  T  ,  responses  u  and  u  were  obtained 

as  shown  in  Fig.  A. 5.   The  displacement  response  u  agrees  well  with 

the  exact  response  shown  in  Fig.  A. 3  for  At  equal  to  0.05  T   and 

0.01  T   but  shows  considerable  error  for  At  equal  to  0.5  T  .   These 

and  similar  results  for  other  values  of  At  indicate  that  the 

approximate  response  u  will  agree  reasonable  well  with  the  exact 

response  provided  At  is  assigned  a  value  less  than  0.2  T  .   The 

velocity  response  u  in  Fig.  A. 5  shows  an  erroneous  oscillation  for 

At  equal  to  0.05  T„  and  0.01  T, .   However,  if  a  smooth  curve  is  drawn 

II 

through  these  oscillations,  it  would  agree  reasonable  well  with  the 
exact  relation  shown  in  Fig.  A. 3.   The  velocity  response  in  Fig.  A. 5 
for  At  equal  to  0.5  T  is  greatly  in  error  which  again  shows  that  At 
must  be  taken  much  smaller  than  0.5  T  . 

Further,  let  us  consider  the  effect  of  the  value  of  n  on 
dynamic  response.   Using  K   equal  to  100,   At  equal  to  0.01  T  ,  and 
n  equal  to  2,  10,  and  40,  responses  u  and  u  were  obtained  as  shown 
in  Fig.  A. 6.   The  approximate  displacement  response  for  all  three  values 
of  n  agrees  very  closely  with  the  exact  response;  however,  the 
approximate  velocity  response  agrees  closely  with  the  exact  response 


A-4 


f 


only  for  the  case  of  n  equal  to  40.   Considerable  disagreement  in  the 
velocity  response  is  noted  for  n  equal  to  10  and  an  even  larger 
disagreement  is  noted  for  n  equal  to  2. 

Finally,  let  us  consider  the  variations  of  particle  velocity 
and  longitudinal  stress  with  position  along  both  bars  for  the  case  of 
n  equal  to  10,   K   equal  to  100  (y  =  1) ,  and  At  equal  to  0.05  T  . 
The  resulting  variations  along  both  bars  for  instantaneous  times  t 
equal  to  0,  0.25  T  ,  0.50  T  ,  0.75  T  ,  1.00  T  ,  and  1.25  T  are  shown 
in  Fig.  A. 7.   The  ordinates  representing  particle  velocity  and  stress 
have  been  normalized  by  dividing  by  the  initial  bar  velocity  (0.1)  and 
the  exact  intensity  of  the  stress  wave  (V   p  E  =  0.2  p  E),  respectively. 
Comparing  these  variations  with  the  exact  solution  shows  that  reasonably 
accurate  results  can  be  obtained  by  the  approximate  method  employing  an 
impact  spring. 

Based  on  the  above  discussion,  the  effects  of  parameters  K  , 
At,   and  n  on  the  dynamic  response  of  two  colliding  rods  can  be 
summarized  as  shown  in  Table  A.l. 
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(a)    APPROACH 


(b)     COLLISION 
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(c)     SEPARATION 


FIG.  A.l   LONGITUDINAL  COLLINEAR COLLISION  OF  TWO  RODS 
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FIG.  A. 3   EXACT  SOLUTION  OF  RELATIVE  DISPLACEMENT  AND  RELATIVE  VELOCITY 
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FIG.  A. 4.1   EFFECT  OF  IMPACT  SPRING  STIFFNESS  ON  RELATIVE  RESPONSE  OF 
RODS  (Kx=  10,  102  and  103) 
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FIG.  A. 4. 2   EFFECT  OF  IMPACT  SPRING  STIFFNESS  ON  RELATIVE  RESPONSE  OF 
RODS  (K1  =  104  and  105) 
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FIG.  A. 5   EFFECT  OF  TIME  INTERVAL  ON  RELATIVE  RESPONSE  OF  RODS 
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FIG.  A. 6   EFFECT  OF  ELEMENT  NUMBER  ON  RELATIVE  RESPONSE  OF  RODS 
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APPENDIX  B 
FREE  VIBRATION  OF  A  SINGLE  SPRING-MASS  SYSTEM  WITH  COULOMB  FRICTION 

Consider  the  rigid  mass  m  shown  in  Fig.  B.l  which  is 

attached  to  a  linear  elastic  spring  of  stiffness  k  and  is  supported  on 

a  surface  which  can  develop  Coulomb  friction  having  a  maximum  absolute 

value  of  Vmg  (V  =  Coulomb  friction  coefficient) .   If  the  mass  is 

displaced  to  the  right,  elongating  the  spring  by  an  amount  x  ,   where 

x  »  A  and  where 
o 

a  -  t  •-  ys. 

,   and  it  is  then  released  from  rest,  the  free  vibration  equation  of 
motion  becomes 

•  •  • 

mx+kx=-F  sign(x)  (B-2) 

Introducing  two  new  variables ,  namely 


x  =  x  +  A 


x2  =  x  -  A 


(B-3) 


,  Eq.  (B-2)  can  be  written  in  the  double  form 


mx+kx=o  x>o 


mx+kx=o  x<o 


(B-4) 


,  having  solutions 


x,  =  A  cos  (0)  t  -  Y)      x  >  o 
1  n 


x„  =  A1  cos  (0)  t  -  Y')     x  <  o 
2  n 


(B-5) 


where  00  =  /k/m  and  A,  A',  y,      and  Y'   are  arbitrary  constants.   The 
n 

solution  of  Eq.  (B-2)  can  now  be  written  as 

B-l 


X  =  -  A  +  A   COS  (d)  t  -  Y)        X  >  o 

n 

(B-6) 
X  =  A  +  A'  COS  (OJ  t  -  Y')         X  <  o 

n 


As  an  example,  let  x   equal  12A.   The  resulting  motion  of  the  mass  as 

represented  by  Eq.  (B.6)  is  that  motion  shown  in  Fig.  B.2.   It  is  well 

known  that  the  frequency  of  free  vibration  of  the  system  with  Coulomb 

friction  is  the  same  as  the  free  vibration  frequency  with  no  damping. 

Further,  it  is  well  known  that  the  amplitude  of  oscillation  diminishes 

by  an  amount  4A  with  each  full  cycle.   The  block  finally  comes  to  rest 

in  the  first  extreme  position  which  is  less  than  A.   For  this  example, 

the   final   rest  is  reached  at  time     t  =  6tt/u)   . 

n 

The  above  example  problem  has  also  been  solved  using  the 
mathematical  model  and  analytical  procedures  defined  in  Chapter  II  and 
III,  respectively.   The  numerical  data  used  in  this  approximate  solution 
were  the  following: 

2 
m  =   80/386.4  lbs  sec  /inch 

k  =  400  lbs/inch 

V  =   0.5 


1    /(400) (386) 

:   =   —   / — =   6.98  c/s  (B-7) 

n     2fT   / 


80 


OJ   =   27T  (6.98)   =   43.8  c/s 
n 


.      F      (0.5)  (80)  . 

A  =  —  =  rr^ =   0.1  inch 

k        400 


x   =   12  A  =   1.2  inch 
o 


In  the  numerical  calculation,  it  is  convenient  to  introduce  the  elastic 
deformation 

B-2 


e      F 
u   =  —  (B-8) 

ke 


e 

where  F  =  Vmg  and  k   is  the  elastic  stiffness  of  the  Coulomb 

friction  elastoplastic  model  defined  in  Fig.  2.4.   The  constant 
acceleration  integration  procedure  with  equilibrium  interation  is  used 
for  the  analysis.   The  tolerances  controlling  the  equilibrium  iteration 
were  specified  as 

<5  R    <_  0.1    ;    <S  R(l)  <_  0.01  (B-9) 

where  6  R   and  6  R     represent  the  residual  forces  at  time  t  and 

during  the  ith  iteration,  respectively. 

e 
First,  let  us  examine  the  effect  of  the  value  of  u   on  res- 
ponse.  Using  a  time  interval  At  equal  to  0.05(2tt/(jo  )   throughout  the 
calculations,  the  displacement  response  shown  in  Fig.  B.3  is  obtained. 

Clearly,  the  approximate  solution  approaches  the  exact  solution  with 

e 
decreasing  values  of  the  ratio  u  /A.   Quite  good  agreement  is  shown 

when  this  ratio  is  assigned  the  value  0.1,  particularly  for  the  higher 

amplitude  oscillations. 

Now  let  us  examine  the  effect  of  the  value  of  At  on  response. 

For  this  examination,  let  u   equal  0.1  A  and  compare  the  resulting 

responses   for     At     equal   to   0.05(2tt/oj   ),      0.10(2tt/w   )/      and     0.166(2tt/w   ) 

n  n  n 

as  shown  in  Fig.  B.4.   It  is  quite  apparent  from  these  results  that  the 
approximate  solution  is  improved  by  decreasing  the  time  interval  At. 
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